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PART 1 


1.1 Draw a block diagram for the control system generated when a human 
being steers an automobile. 



1.2 From the given figure specify the devices 



































Solution: 

Loop 110 







Sequential indicating controller. Y represents mathematical symbols. 


Emits signal to valve which controls the flow rate of B. 


-t>o- 


- solenoidal valve. It operates on electric signal. 


Loop 103 



h- 


Orifice meter. Measures flow rate. 



Flow transmitter.lt converts flow to pneumatic signal of 20 to 100 
Mpa. 


Mr - 


Removes the square root from signal and makes it linear. Input 
pneumatic and output is electrical. 


is 


Loop 104 



Flow transmitter converts flow of A to a pneumatic signal of 20-100 
Kpa. 



Pressure controller. It receives two electrical signals one from flow 

Transmitter of loops 103 and 104.Maintain ratio of A and B in 
reaction chamber. 



IP converter, it converts standard electric signal to pneumatic signal. 




A control valve. Input is a standard pneumatic signal and controls 


flrwi/ rofo 












Loop 101 



Pressure transmitter. Measures pressure and given electric signal of 4 
20 Ma. 



I 


-o- 





Computer. Sends digital data to pressure controller. It is used to set 
the set point and placed in control room. 

Pressures controller - has two inputs one digital signal from 
computer and one electrical signal from pressure transmitter. 

Motor controlled valve,controls the flowrate of C in the chamber. 


Loop 102 


Temperature safety. 

m- 

Sequential indicating controller. 

Loop 112 


Sequential indicating controller 

I><1 

Solenoidal valve. Controls flow rate of vent gases 


Inversion by partial fractions: 
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5 2 X(5) - sx(O) -x (0) + s X(5) -x(0) + 2 /( 5 ) = — 
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= (s 2 + 5 + l)X(5) =- 
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X(s)= T^- n 

s(s +5 + 1) 

Now, applying partial fractions splitting, we get 
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dr 2 dx 

b) —5- + 2 — + x = 1 jc(0) = jc (0) = 0 

dt dt 


when the initial conditions are zero, the transformed equation is 


(5 2 + 5 + l)X(5) =- 
5 
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X(s) 


.S'(.S' 2 +5 + 1) 


1 A Bs + C 

- 9 - — - 1 - 9 - 

5(5" +5 + 1) 5 5" +25 + 1 

1 = 4++2s+1)+£s 2 +Cs 

0 = A + B(by equating the co - effecient of s 2 ) 
0 = 2 A + C(by equating the co - efficients of 5 ) 


1 = A (by equating the co - effecients of const ) 
A + B = 0 
B = -1 
C = -2 A 

A- l,B - -1,C - -2 


X(s) = -— 2 S + 2 

5 5 2 + 25 + 1 

r 1 m5)} = r 1 (l- ( ; + 1 |; 1 

[5 (s+l) , 


{x(0} = i-r 1 j^- + r ^— 1 

[5 + 1 (5 + 1)- J 
{X(t)} = l-e\\ + t) 

dx ^ dx * 

3.1 C —-r- + 3—- + x = 1 x(0) = x(0) = 0 

dr dt 


by Applying laplace transforms, we get 

= (5 2 + 35 + 1)X(5) = - 
5 


X(5) = 


5 ( 5 ^ +35 + 1) 



X(s) = - + 
s 


Bs + C 
s + 3s +1 


1 = A(s 2 + 3s + 1) + Bs 2 + Cs 
0 = A + B(by equating the co - effecient of s 2 ) 

0 = 3 A + C(bv equating the co - effecients of s) 


1 = A(by equating the co - effecients of const) 
A + B = 0 
B = -1 

C = -3A = -3 
A = \,B = -1,C = -3 


Z-'{X(5)} = L- 1 |i 


5 + 3 | 

5 2 +35 + lj 


Z-‘{X(5)} = Z- 1 
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X(t) - I - e 2 (Cos 2^ + -^=sinh 1 

2 S 2 


3.2(a) 

dr 4 d 3 x 

^ + ffb = C ost;x(0) = x (0) = x (0) = 0 
dt dt 


Jt“(0) = l 



Applying Laplace transforms, we get 


s 4 X(s) - s 3 x( 0 ) - sV( 0 ) - ,vx ( 0 ) - X ( 0 ) + s 3 X(s) - s 2 x( 0 ) - sjc’(O) - jc"(0) = — r- — 

.S' - +1 


X(s)(s 4 +s 3 )-(s + 1) = -^- 

s +1 


X(5) 


r 
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+ l) + 0 + l) 


/s 4 +s 2 
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S + S 3 +S+S 2 +l 

sV+l)(s + l) 


s s 4" 2 s +1 
s 3 (s 2 +1)(J + 1) 


,v + s + 2 s + 1 A B C D Es + F 

s\s 2 +1)(5 + 1) 5 s 2 S 3 S + \ 5 2 +l 


s 3 + s 2 + 25 + 1 = As 2 (s + 1)(V + 1) + Bs(s +l)(s 2 +1) + c(s +1)(> 2 +1) + Ds\s 2 +1) + (Es + F> 3 (s +1) 


A+B+E=0 equating the co-efficient of s 5 . 
A+B+E+F=0 equating the co-efficient of s 4 . 
A+B+C+D+F=0 equating the co-efficient of s 3 . 
A+B+C=0 equating the co-efficient of s 2 . 
B+C=2 equating the co-efficient of s. 
A+B+E=0 equating the co-efficient of s 2 . 

C=1 equating the co-efficient constant. 


C=1 

-B=-C+2=l 

A=l-B-C=-1 

D+F=0 

E+F=0D+E=1 

D-E=0 

2D=1 


A=-l; B=l; C=1 
D=l/2; E=l/2; F =-1/2 



L~' {(*)} = L~' 


fi>4 + *3 + M 

[ S S S 5 + 1 s 2 +1 J 



1/2(5-!) ] 

5 2 +1 J 


{X(t)} = -l + t + — + — e ‘ + — Cost - — Sint 
2 2 2 2 


^f + ^W+2 1 q{0) = A-q\0) = -2 
at at 

applying laplace transfonns,we get 

s 2 Q(s) - sq( 0) - q'( 0) + sQ((s) - g(0) = \ \ 

s s 

2 (1 A 

Q(s)(s 2 +s)-4s+ 2-4 = — - + 1 

5 U ) 




3.3 a) 


3s \ 1 


(5 2 +l)(r+4) 3ls z +l s 2 + 4 


,2 , ! 2 


s 2 + l 2 s 2 + 2 2 


s 2 +l 2 s 2 + 2 2 


= Cost - Coslt 


A B + C 


2 , o2l _ + „2 


5(5 2 -25 + 5) 5(5-ir+2 2 5 5 2 -25 + 5 


A+B=0 

-2A+C=0 

5A=1 

A=l/5 ;B=-l/5;C=2/5 


We get 


1 1 2-5 

X(s) = -+ —- 

5 5 5-25 + 5 


Inverting,we get 


= — 1 + — e'Sinlt - e'Coslt 
5 2 


1 (1 

= - l + e f —Sinlt - Cos2t 
5 v 2 


3s 2 -s 2 -3s + 2 ABC D 

Cl -T-T- _ -1- T ^ -1-o 



As(s - 1) 2 + B(s - 1) 2 + Cs\s - 1 ) + Ds 2 = 3s 2 - s 2 - 2s + 2 


A(s 3 -2 s + s) + B(s 2 -2^ + 1) + C(s 3 - 5 2 ) + Ds 2 = 3s 3 - s 2 - 3s + 2 


A+C=3 

-2A+B-C+D=-1 

A-2B=-3 

B=2; 

A=2(2)-3=l 

C=3-l=2 

D=2(l)-2+2-l=l 


We get X(s) = - + 4 + ^T+ , 

S S 5 + 1 (5-1) 

By inverse L.T 

L~ l [X(t)\=l + 2t + 2e l +te‘ 
r 1 [X(t)] = l + 2t + e , (2 + t) 


3.4 Expand the following function by partial fraction expansion. Do not evaluate 
co-efficient or invert expressions 


X(5) 


_ 2 _ 

(5 + l)(5 2 +l) 2 (5 + 3) 


X(5) 


A Bs + C Ds + E F 

-1-^-1-^- 7 -l- 

5 + 1 5" + 1 (5" + l) - 5 + 3 


= A(s 2 + 1) 2 (5 + 3) + (Bs + C)(5 +1)(5 + 3)(5 2 + 1) + (Ds + E)(5 + 1)(5 + 3) + F(s + 1)(5 2 + 1) 2 


= A(s 4 + 25 2 + 1)(5 + 3 ) + (Bs + c )( 5 2 + 45 + 3 )( 5 2 + 1 ) + (Ds + E)(s 2 + 45 + 3 ) + F(s + 1 )( 5 2 + 45 + 1 ) 



= 5 s (ri + 5 + F) + s 4 (3A + C + 4B + F) + s 3 (2A + B + AC + 3B) + s 2 (6A + C + 4B + 3C) + s(A + 4C + 35 
+ 4E + F) + 3A + 3AC + 3E + F = 2 
A+B+F=0 

-3A+C+4B+F=0 
2 A+B+4C+3 B=0 
6 A+C+4B+3 C=0 
A+4C+3B+3D+4E+F=0 
3A+3C+3E+F=2 

by solving above 6 equations, we can get the values of A,B,C,D,E and 

X(s)=^ --- T . 

5 3 (5 + l)(5 +1) (s + 3) 3 


V .._ABC D E F 

A(s) —-1-—4- r- 4-1-1-h 

5 S 2 5 3 5 + 1 5 + 2 5 + 3 


by comparing powers of s we can evaluate 
c)X(5) =--- 

5(5 + 2)(5 +3) (5 + 4) 


G H 

(5 + 3) 2 + (5 + 3) 3 

A,B,C,D,E,F,G and H. 


w , A B C D 

X(s) =-+-+-+- 

5+1 5+2 5+3 5+4 

by comparing powers of s we can evaluate 


A,B,C,D 


3.5 a) X{s) 
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B/2=l/2 *-3/2=-l; 
B=-2; 

C= -3/2+2=l/2 


w , _ 1 2 If 1 > 

X(s) — -1— - 

s s +1 2^0.55 + 1 


= L 1 (X (s)) = x(t) = 1 - 2 e 1 + e 


-it 


dx 

b) — + 2x = 2; x(0) = 0 
dt 


Applying laplace trafsorms 


sX(s) - x(0) + 2X (s) = 2/s 


L\X(s)) = 


s(s + 2) 


L\X(s)) = 2L- 


s(s + 2) 


=L~\X{s)) = 2L~ 


1/2 1/2 


s s + 2 


= 1 - e 


-it 


3.6 a) 7 ( 5 ) = 


5 + 1 


s + 2 s + 5 



5 + 1 

5 ~ + 2 s + 5 


= 7(5) = 


5 + 1 

(5 +1) 2 +4 


= L- l (Y(s)) 



5 + 1 

(5 +1) 2 +4 


using the table,we get 

Y(t ) = e'Coslt 


5 2 + 25 
b) 7(5) = ^ 

5 

7 ( 5 ) = 4+4 

5 5 

Y(t)= r'(7(5)) = t + t 2 


c) 7(5) 


25 

(.S'-l ) 3 


25-2 + 2 

(.S'-l ) 3 


2 2 

(.S'-I) 2 + (.S'-1) 3 


Y(t)=L l 


(s-\y 


+ r l 


(5-iy 


= 2 (te +—e r = e (f + 2 1) 



3.7a) Y(s) 


1 _ As+ B Cs + D 

(r+1) _ (r+l) + (r+l) 


thus (As + B ) + (Cs + D)(s 2 + 1) = 1 
= Cs 3 +Ds 2 +(A + C)s + (B + D) = 1 


C=0,D=0 
Also A=0;B=1 

w , 1 1 A B C D 

(s 2 +l ) 2 (s + i) 2 (s-i ) 2 (s + i ) 0 + if ( s-i) ( s-if 


A(s + i)(s - i ) 2 + B(s - i ) 2 + C(s - i)(s + i ) 2 + D(s + i ) 2 = 1 


(A + C)s 3 + (-Ai + B + Ci + D)s 2 +(A-2 Bi + C + 2 Di) + (-Ai -B + Ci-D) = 1 


Thus,A+C=0 

-Ai+B+Ci+2Di=0 ; B=D 
A-2Bi+C+2Di=0 

-Ai-B+Ci-D=l Also D=-Ci;B=-Ci, 

A=-C,C=-i/4 


A=i/4 ; B=-l/4; D=-l/4 


m=——+ 


i /4 - 1/4 -i/4 -1/4 


■ +-+ - 


(5 + i) (s + i ) (s - i ) (s - i ) 
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7(0 = --+-- +-+ ■ 


(s + i) (s + i ) 2 (s-i) (s-i ) 2 


-1/4 

(s-i) 2 



7(0 = i / Ae~ u - 1 / Ae ~ u -\/4e' Y -1 /4fe' Y 


7(0 = 1 / 4(ie-" - te~ u - ie" - te“) 

7(0 = 1 / 4(/(Co.st - /.S7« 0 - t(Cost - i Sin t ) - i(Cos t + iSin t) - t(Cos t + i Sin t )) 
7(0 = 1 / 4(2 Sin t-2t Cos t) 

Y ( t ) = 1 / 2 ( t - t Cos t ) 


3.8 f(s) 
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, A B 

~ f( s )-— + — + 

s s 


C 

5 + 1 


= ,4(5 + l) + 55(5 + l) + Cr =1 


Let s=0 ; A=1 

s=l; 2A+B+C=1 
s=-l: C=1 

B=-l 


„ , _ 1 1 

f( s ) - - + - + 

s s 


1 
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PROPERTIES OF TRANSFORMS 


4.1 If a forcing function f(t) has the laplace transforms 


m = -+ 


1 ^-2.s -3s 

1 e -e e 


i 3 a a -25 

1-e e -e 

-+-r- 



m = L~ x {/($)} = [u(t) - u(t - 3)] + [(t - 1 )u(t -1) - (f- 2)u(t - 2)] 


- u(t) + (t - 1) u{t - 1) - (t - 2)u(t - 2) - u(t - 3) 
graph the function f(t) 


a j 

/ . 


\ 
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f i6io9€-- 


4.2 Solve the following equation for y(t): 

\y(r)dr = d ^~ y(0 ) = 1 

Taking Laplace transforms on both sides 

L \ y{Tut)=L W\ 

-.y(s)=s.y(s)-y(0) 

s 










-■y(s) = s.y(s)-l 
s 



y(t)=L l {v(s)}= V' !^| cosh(t) 

4.3 Express the function given in figure given below the t - domain and the 
s- 

domain 


This graph can be expressed as 

= {u(t - 1 ) — u(t — 5)} + {(t- 2)u(t — 2) — (t- 3 )u(t - 3)} + {u(t — 5) - (t — 5)u(t - 5) + (t - 6 )u(t - 6)} 
f(t)= u(t - 1) + (t - 2)u(t -2) — (t- 2)u(t -3 ) — (t- 5)u(t — 5) + (t — 6 )u(t - 6) 


f(s)=L{f(t)} = — + ^ 
s s 


-6 s 


- + - 
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4.4 Sketch the following functions: 
/(7) = w(7) - 2u(t -1) + u(t - 3) 


Z - 



v-4 -4 C*} 


f (t) = 3tu(t) - 3 u(t -1 ) + u{t - 2 ) 





















V \ fc\ u (? ^ 2+- k? 


4.5 The function f(t) has the Laplace transform 

f(S)= (1 - 2e~ s +e~ 2s )/ s 2 


obtain the function f(t) 


f(s) = 


1 - 2e~ s + e~ 


and graph f(t) 


-2s 


f{t)-L~ l {f(s)} = - (t -1 )u(t -1) + tu(t) - {{t -1 )u(t - 1>- (t - 2)u(t - 2)] 


= tu(t) - 2(t -1 )u(t -1 ) + (t - 2)u(t - 2) 

4.6 Determine f(t) at t = 1.5 and at t = 3 for following function: 


/( t ) = 0.5«(Y) - 0.5 u(t -1) + (t - 3 )u{t - 2) 









At t = 1.5 


/( t ) =0.5 u(t) - 0.5 u(t -1 ) + (t - 3 )u(t - 2) 
/(1.5) =0.5 u(t) - 0.5 u(t - 1) 

/ (1.5) =0.5-0.5 = 0 
At t = 3 

/(3) =0.5—0.5+(3—3)= 0 


RESPONSE OF A FIRST ORDER SYSTEMS 


5.1 A thermometer having a time constant of 0.2 min is placed in a 
temperature bath and after the thermometer comes to equilibrium with 
the bath, the temperature of the bath is increased linearly with time at the 
rate of I deg C / min what is the difference between the indicated 
temperature and bath temperature 

(a) 0.1 min 

(b) 10. min 

after the change in temperature begins. 

© what is the maximum deviation between the indicated temperaturew 
and bath temperature and when does it occurs. 

(d) plot the forcing function and the response on the same graph. After the 
long enough time buy how many minutes does the response lag the input. 


Consider thermometer to be in equilibrium with temperature bath at 
temperature Xs 
X(t) = X s +(l°/m)t,t>0 

as it is given that the temperture varies linearly 
X(t)-Xs = t 


Let X(t) = X(t) - Xs = t 



Y(s) 

Y(s) 


= G(s).X(s) 

11 ABC 

— - y — - 1 - 1 - z 

l + ws 1 + w s s 


A = T 2 B = -t C= 1 


^ 1 2 
1 + TS S S 

Y(t ) = re~' /r - r + t 

(a) the difference between the indicated temperature and bath temperature 
at t = 0.1 min = X(0.1)_ Y(0.1) 

= 0.1 - (0.2e' 01/a2 - 0.2+0.1) since T = 0.2 given 
= 0.0787 deg C 

(b) t = 1.0 min 

X(l) - Y(l) = 1- (0.2e‘ 1/0 '2 - 0.2 +1) = 0.1986 

(c) Deviation D = -Y(t) +X(t) 

= -xe" t/T +T =x (-e’^+l) 

For maximum value dD/dT = x (-e‘ t/T +(_-l/T) = 0 
-e _t/ = 0 

as t tend to infinitive 
D = x (-e‘ t/T +(_-l/T) = x =0.2 deg C 


5.2 A mercury thermometer bulb in Vi in. long by 1/8 in diameter. The 
glass envelope is very thin. Calculate the time constant in water flowing 
at 10 ft / sec at a temperature of 100 deg F. In your solution, give a 
summary which includes 


(a) Assumptions used. 

(b) Source of data 

(c) Results 



T = mCp/hA 


(pAL)C p 
h(A + kDL) 


Calculation of 

NU d ~ = CR';m n 

_ Dvp _ (1/8 *2.54 * 10~ 2 )(10*0.3048)10 3 _ n£nn A 

IvU ; - - - - 7U / / . I 

d jU 10' 3 

Pr = —^ = 4.2 KJ/KgK 
K 

Source data: Recently, Z hukauskas has given c,m ,q,n values. 
For Re = 967704 
C = 0.26 & m = 0.6 

Nu d = hD/K =0.193 (9677.4)*(6.774X10' 3 ) = 130 
.h = 25380 


5.3 Given a system with the transfer function Y(s)/X(s) = (Tis+1)/(T2S+1). 
Find Y(t) if X(t) is a unit step function. If T 1 /T 2 = s. Sktech Y(t) Versus 
t/T 2 . Show the numerical values of minimum, maximum and ultimate values 
that may occur during the transient. Check these using the initial value 
and final value theorems of chapter 4. 


7(5) 


?> + l 

X 2 S + 1 


X(s) =unit step function = 1 X(s) = 1/s 



Y(s) 


T,s +1 A B 

— 1 -= — +- 

s(T 2 s + 1 ) s iT 2 s 


A = 
Y(s) 


IB = Ti - T 2 

1 T - T 
1 1 1 2 

s i + r 2 s 


7(0= 1 + 


^1 T 2 „-t/T 2 


If T 1 /T 2 = s then 
7(0= l + 4e~‘ IT2 


Let t/T 2 = x then 7(0 = l + 4e ' 

Using the initial value theorem and final value theorem 
LimY(T) = LimsY(s) 

r -> 0 s -*- 00 


r+ A 

T r + I U T 

= Lim — -= Lim -A 

S —>00 T 2 s +1 S'—> °° 1 

^2 ' 


7 


Lim Y (T) = Lim sY (.v) = Lim ^ + ^ = 1 

S—> 00 


r->o 


J 2 5 + 1 


5 


Figure: 




P'Tlobl^'K 5- 3 


5.4 A thermometer having first order dynamics with a time constant of 1 
min is placed in a temperature bath at 100 deg F. After the thermometer 
reaches steady state, it is suddenly placed in bath at 100 deg F at t = 0 and 
left there for 1 min after which it is immediately returned to the bath at 
100 deg F. 

(a) draw a sketch showing the variation of the thermometer reading with 
time. 

(b) calculate the thermometer reading at t = 0.5 min and at t = 2.0 min 


X(s) s + l 


1 min) 


(S) = 10 


5 


5 














7(s) = 10 


1 - e~ s 


1 


7(s) = 10 

s(s + 1) s(s + 1) 
7(0 = 10(1-0 t<\ 


7(0 = 10( (1 - e~‘) - (1 - e" (M) ) )t > 1 
At t = 0.5 T = 103.93 
At = 2 T =102.325 


p7l© 








5.5 Repeat problem 5.4 if the thermometer is in 110 deg F for only 10 sec. 

If thermometer is in 110 deg F bath for only 10 sec 
T = 110 - 10O 60 












0 < t < 10 sec & T = 60 sec 
T(t = 10sec) = 101.535 
T =100 +1.535e _( '“ 10)/60 t>10sec 
T(t=30sec) = 101.099 deg F 

T(t=120sec) = 100.245 deg F 

S'-S * 



/ . 

5.6 A mercury thermometer which has been on a table for some time,is 
registering the room temperature ,758 deg F. Suddenly, it is placed in a 400 
deg F oil bath. The following data are obtained for response of the 
thermometer 


Time (sec) 

Temperature, Deg F 

0 

75 

1 

107 

2.5 

140 

5 

205 

8 

244 

10 

282 

15 

328 

30 

385 


Give two independent estimates of the thermometer time constant. 


















T = 


t 

325 


400 -T 


From the data, average of 9.647,11.2,9.788,10.9,9.87,9.95, and 9.75 is 10.16 sec. 


5.7 Rewrite the sinusoidal response of first order system (eq 5.24) in 
terms of a cosine wave. Re express the forcing function equation (eq 5.19) 
as a cosine wave and compute the phase difference between input and 
output cosine waves. 


r(g)=^-r(g)= 2 Ao) 2 1 

zs +1 s + co , 1 
s + - 


splitting into partial fractions then converting to laplace transforms 


„. . Aacor A . . .. 

^(0= -r“ 2 —~ e + i = „ sin (at + </>) 


t 2 co 2 +1 


rV+1 


where (f» = tan" (cox) 


As t —>oc 


^(0[s = , A 1 sin {cot + (t)) = — A cos ) 

Jr~co~ +1 Vr"®‘+1 V 2 J 


Y(t) = Asi\\{cot + 0) = Acos - cot 

v 2 


Y(t) = A cos cot - 

V 2 


The phase difference = </)- — - 

2 V 2 ) 



5.8 The mercury thermometer of problem 5.6 is allowed to come to 
equilibrium in the room temp at 75 deg F.Then it is immersed in a oil 
bath for a length of time less than 1 sec and quickly removed from the 
bath and re exposed to 75 deg F ambient condition. It may be estimated 
that the heat transfer coefficient to the thermometer in air is 1/ 5 th that in 
oil bath.If 10 sec after the thermometer is removed from the bath it reads 
98 Deg F. Estimate the length of time that the thermometer was in the bath. 

t < 1 sec T x = 400 - 325e-' l/r 

Next it is removed and kept in 75 Deg F atmosphere 

Heat transfer co-efficient in air = 1/5 heat transfer co-efficient in oil 


hair 1/5 h 0 ii 


mC 

t= - T ojl = 10 sec 


hA 


T air = 50 sec 


T f = 75 + {T x -75)e"' 2/50 
T f = Final Temp = 98 deg C 
98= 75 + (325 -325A' l/10 )A 10/50 

e~ ,no =0.91356 
ti = 0.904 sec. 


5.9 A thermometer having a time constant of 1 min is initially at 50 deg C. 
it is immersed in a bath maintained at 100 deg C at t = 0 . Determine the 
temperature reading at 1.2 min. 

x = 1 min for a thermometer initially at 50 deg C. 

Next it is immersed in bath maintained at 100 deg C at t = 0 


At t = 1.2 



7(f) =A(l - e~ t,T ) 


7(1.2) = 50(1 -e~ 12/1 ) + 50 
Y(1.2) = 84.9 deg C 

5.10 In Problem No 5.9 if at, t = 1.5 min thermometer having a time 
constant of 1 minute is initially at 50 deg C.It is immersed in a bath 
maintained at 100 deg C at t = O.Determine the temperature reading at t = 
1.2 min. 

At t = 1.5 

7(1.5) = 88.843°C 

Max temperature indicated = 88.843 deg C 
AT t = 20 min 
T = 88.843 - 13.843(1-e" 18 ' 8/1 ) 


T = 75 Deg C. 


5.11 A process of unknown transfer function is subjected to a unit impulse 
input. The output of the process is measured accurately and is found to be 
represented by the function Y(t) = t e *. Determine the unit step response in 
this process. 

X(s) = 1 Y(t) = te' 1 




i 

(.v + lf 


G(s ) 


TOO 

X{s) 


1 

(s + 1) 2 


For determining unit step response 


Y(s) 


1 

(s + 1) 2 



1 


c 


Y(s) 


1 _ A B _ 

(.S’ +1) 2 s 5 + 1 (5 +1) 2 


A = 1 B = -1 C = -1 


Y(s) = 


]_ _ 1 1 

5 5 + 1 (5 + l) 2 


7(0= 1 -e~‘ -te~‘ 


Response of first order system in series 


7.1 Determine the transfer function H(s)/Q(s) for the liquid level shown in 
figure P7-7. Resistance Ri and R 2 are linear. The flow rate from tank 3 is 
maintained constant at b by means of a pump ; the flow rate from tank 
3 is independent of head h. The tanks are non interacting. 



Solution : 

• \ • ¥ • > t VX- 1 Uf I i W/^'ML / % 



PlGiO^F P7-1 


A balance on tank 1 gives 


. dh , 

q ~ q ' =A >^~ 
dt 


where hi = height of the liquid level in tank 1 


similarly balance on the tank 2 gives 

. dh., 

d\ ~ di = ^2 


dt 


and balance on tank 3 gives 






















dl do 


A 


dh 

dt 


here q x 


A_ 

R, 


di 



b 


So we get 


d - 


K_ 

R, 


A 


dh l 

dt 


hi_hi_ 

A R 2 


a 2 


dh 2 

dt 


h 2 

A 


-b 


A 


dh 

dt 


writing the steady state equation 


ds 


h 


Is 

A 


A 


dh is 

dt 


0 


As 

A 


h 

n 2S 

r 2 


A 


dh~, s 

dt 


R 2 

Subtracting and writing in terms of deviation 


Q- 


H_ 

A 


A 


dH x 


dt 



dH 2 

dt 


i?i r 2 

H 2 _ dH 
R 2 dt 

where Q = q -qs 
Hi= h r hi S 

Hl= h2-ll2S 

H = h - h s 

Taking Laplace transforms 


Q(s) - 


H x (s) 

*1 


A x s H x (s) 


H,{s) H 2 (s) 

R \ R 2 


A 2 s H 2 (s) 


( 1 ) 


( 2 ) 


= sH(s) .(3) 

R 2 

We have three equations and 4 unknowns(Q(s),H(s),Hi(s) and H 2 (s). So we 
can express one in terms of other. 

From (3) 


H 2 (s) 

R x A 2 s 


H 2 (s ) 


(4) 







H 2 (s) 


Rl - Hl ^ where r^IiA 
Ri(T 2 s + 1 ) " - 


(5) 


From (1) 


H { (s) 


R,Q(s) 

{t x s + 1 ) ’ 


T i = Ri A -( 6 ) 


Combining equation 4,5,6 


H (s) 


Q(s) 

(A 3 s)(r t s +1 )(t 2 s + 1) 


H(s) 1 

Q(s) (AjsX^s + 1)(t 2 s + 1 ) 


Above equation can be written as 

i.e, if non interacting first order system are there in series then there overall 
transfer function is equal to the product of the individual transfer function in 
series. 


7.2 The mercury thermometer in chapter 5 was considered to have all its 
resistance in the convective film surrounding the bulb and all its 
capacitance in the mercury. A more detailed analysis would consider both 
the convective resistance surrounding the bulb and that between the bulb 
and mercury. In addition , the capacitance of the glass bulb would be 
included. 

Let 

Ai = inside area of bulb for heat transfer to mercury. 

Ao = outside area of bulb, for heat transfer from surrounding fluid. 

.m = mass of the mercury in bulb. 
m b = mass of glass bulb. 

C = heat capacitance of mercury. 





Cb = heat capacity of glass bulb. 

.hi = convective co-efficient between the bulb and the surrounding fluid. 

.ho = convective co-efficient between bulb and surrounding fluid. 

T = temperature of mercury. 

Tb = temperature of glass bulb. 

Tf = temperature of surrounding fluid. 

Determine the transfer function between Tf and T. what is the effect of bulb 
resistance and capacitance on the thermometer response? Note that the 
inclusion of the bulb results in a pair of interacting systems, which give an 
overall transfer function different from that of Eq (7.24) 


Writing the energy balance for change in term of a bulb and mercury 




Input - output = accumulation 


hMTf-T^MTh-T) = m b C b ^ 


dT 

h t A i {T b -T)-Q = m C — 

dt 


Writing the steady state equation 


Mo (T fs -TJ~hA(T bs ~T) 


m h C b 


dt 


0 


KMT* -T s ) = 0 







Where subscript s denoted values at steady subtracting and writing these 
equations in terms of deviation variables. 


hMT f -T t )-lhMT„-T m ) = 

dt 

h,A,(T t -T.>0= mC 

dt 

Here Tp = Tf - Tfs 
Tb = Tb - Tbs 
Tm= T - T s 

Taking laplace transforms 
hMT^-T^yh.MTB-TJ = m b C b T B (s) —(1) 
And h i A i (T B (s)-T m (s)) = mC sT B (s) .(2) 


= h 0 A 0 (T F (s)-T B (s))-mCST m (s) = m h C b sT B (s ) 


From (2) we get 


T b (s) = T m {s) 


mC 

~h~A 


5 + 1 


= T J S ) (V + 1 ) 


Where r, = 


mC 

h,A. 


Putting it into (1) 


T f {s) - T m (s) 


TYlC 

(r ; .5 + l))(r 0 5 + l) + --—5 
rl 0 A 0 


= 0 




= TJs) = TJs) 


( r,s + l))(r 0 s +1) + 


mC 

h 0 A 0 


TJs) 

TJs) __„2 


1 


z , mC 

TiT 0 S +(T i +T 0 +——)s + l 


A 0 


TJs) _ 1 _ 

t.t 0 s 2 + (r. + r 0 + )s +1 

0 0 h A 

n Q /i Q 


Or we can write 


TJs) 

TJs) 


_1 

T t T 0 S 2 + (T, +T 0 + 


mC 

-> + l 


mC , 

t. = - and r n 

hj 


m h c b 

h 0 A 0 


We see that a loading term mC/ hoAo is appearing in the transfer function. 
The bulb resistance and capacitance is appear in r 0 and it increases the 
delay i.e Transfer lag and response is slow down. 


7.3 There are N storage tank of volume V Arranged so that when 
water is fed into the first tank into the second tank and so on. Each tank 
initially contains component A at some concentration Co and is equipped 
with a perfect stirrer. A time zero, a stream of zero concentration is 
fed into the first tank at volumetric rate q. Find the resulting 
concentration in each tank as a function of time. 


Solution: 




r r v dC ' 

qc t -1 - qC t = v— 

at 

qC,_ X)s - qC is = 0 


c - c 



T = — 

v q) 


V dC t 
q dt 


Taking lapalce transformation 

C (M) (J) - C.(.9) = zsCi(s) 
C (M) (j) = (1 + Ts)Ci(s) 


C,(s) 1 

C t _ j ( 5 ) 1 + r 5 

Similarly 

C^ = C 1 (^ x Q(5) x _ C m (5) ;; C,.(5) 1 

Co(s ) C 0 O) C, (.S’) C,._ 2 (s) C,0) (1 + tj)i 
























Or 


Q(.v) = 1 

Co(s ) (l + rs) w 



CA°) 



T T 

(I+r7F _ (1 + rs) A,_1 





C N (t) =-c 0 





7.4 (a) Find the transfer functions H 2 /Q and H 3 /Q for the three tank system 
shown in Fig P7-4 where Hi,H 3 an d Q are deviation variables. Tank 1 and 
Tank 2 are interacting. 

7.4(b) For a unit step change in q (i.e Q = 1/s); determine H 3 (0), H 3 (oo) 
and sketch H 3 (t) vs t. 

Solution : 

Writing heat balance equation for tank 1 and tank 2 






dC*r 





<7l ^2 


a 2 


dh 2 

dt 


?i 


h x - h 2 


<h 


h-2 

Y 2 


Writing the steady state equation 

(is - Vis = 0 

<h, - <h, = 0 


Writing the equations in terms of deviation variables 

dH x 


Q-Q, = A 


dt 


Q-Q^A^ 


dt 











Taking laplace transforms 

Q(s)-Q 1 (s)= A x sH x {s) 
a (s)-Q 2 (s)= A ]S H 2 (s) 

R t Q x (s)= H 1 (s)-H 2 (s) 

R 2 Q 2 {s)= H 2 (s ) 

Solving the above equations we get 


H 2 {s) _ R, _ 

00) [r,r 2 s 2 + (r, + r 2 + A X R 2 )s + 1] 

Here T x = R ] A i 

Now writing the balance for third tank 


di q 3 


A 


dh 3 

dt 


Steady state equation 


d2S diS ~ ^ 


d 3 = 


h 3 

j; 


02 - 


^3 

R , 


-= A, 


dh 3 

dt 



Taking laplace transforms 


Q 2 (s)-^-=A 3 sH(s) 

02 0) = -j^( T r s + l ) where h= R A 

From equation 1,2,3,4 and 5 we got 

QM __l_ 

00 ) [t x t 2 s 2 + (r, +r 2 +^i? 2 > + l] 

Putting it in equation 6 

H 3 (S) __ i ? 3 

00) [r,r 2 s 2 + (r, + r 2 + A x R 2 )s + l](r 3 s +1) 

Putting the numerical values of Ri R 2 and R 3 and Ai,A 2 ,A 3 

# 3 P) 4 

0p) [4s 2 + 6 s + l](2s + l) 

H 2 (s) 2 

0p) [4s 2 +65 + 1] 

Solution (b) 

0P)=- 

s 



H 3 (s) 


1 4 

5 [4s 2 +65 + l](25 + l) 

From initial value theorem 


H 3 (0) = Lim sH 3 ( 5 ) 

S->00 


Lim --- 

S-> 00 (: 2s + l)(4s 2 +6^ + 1) 


Lim 

S-> 00 


(25 + l)s 3 (4+ + V 2 ) 
s / s 


H 3 (0)= 0 


From final value theorem 

H 3 (co) = Lim sH 3 ( 5 ) 

= Lim -- 

«-*> (25 + l)(45 2 +65 + I) 

h 3 ( 00 ) = 4 









7.5 Three identical tanks are operated in series in a non-interacting fashion 
as shown in fig P7.5 . For each tank R=l, x = 1. If the deviation in 
flow rate to the first tank in an impulse function of magnitude 2, 
determine 

(a) an expression for H(s) where H is the deviation in level in the third 
tank. 

(b) sketch the response H(t) 

(c) obtain an expression for H(t) 

solution : 



writing energy balance equation for all tanks 


q -q x 


A 


dh x 

dt 


<l\ ~<h 


A 


dh 2 

dt 










A 


dh 

dt 


<h 


R 


q 2 



h_ 

R 


So we get 

ds — dis ~ ® 
dis ~dis = 0 
dlS — dis = 0 

writing in terms of deviation variables and taking laplace transforms 

q ( s) _^M = A s H,(s) 

K 


Qi(s) H 2 (s) 
R R 

H 2 (s) H ( j ) 
R R 


A s H 2 (s) 
A s H (s) 


solving we get 


H (s) _ R _ 1 

Q(s) ~(TS + 1) 3 “ (s + lf 

(r 5 +1) 3 (s + 1) 3 



H{t) = L- x {H{ S )} = 2 t -e~ t 


H(t ) = fe 


dt 

= 2t =t 2 

at t = 2 max will occur. 

7.6 In the two- tank mixing process shown in fig P7.6, x varies from 0 lb 
salt/ft 3 to 1 lb salt/ft 3 according to step function. At what time does the 
salt concentration in tank 2 reach 0.6 lb/ ft 3 ? The hold up volume of each 
tank is 6 


( 



Solution 


Writing heat balance equation for tank 1 and tank 2 







Vis -<lys = 0 

Vys -Qcs = 0 


writing in terms of deviation variables 

X(s) - Y(s ) = —sY(s) 

q 

7 ( 5 ) _ 1 _ 1 _V_ 

X(5) (V TS +1 <7 

U ) 

c(s ns) _ X{s) 

K (r 5 +1) (rs+l) 2 

Qs) 1 

X(5) (r 5 +1) 2 



taking laplace transforms 



C(s) = 


(1/4) 

s(s +-) 2 
2 


C(s) = 


f i > 




1.2 


(S +2 ) 


C(s) = 


rn 


1 

U J 


1 

s 

Oil 

2 

2 j 

r n 

s + — 

V 2j 


f ii, 1 

C(t)= \- — te 2 -e 2 


C (0 = 0.61 lb salt / ft 3 


t = 4.04 min 


7.7 Starting from first principles, derive the transfer functions Hi(s)/Q(s) 
and H 2 (s)/Q(s) for the liquid level system shown in figure P7.7. The 
resistance are linear and Ri= R 2 = 1. Note that two streams are flowing 
from tank 1, one of which flows into tank 2. You are expected to give 
numerical values of the parameters and in the transfer functions and to 
show clearly how you derived the transfer functions. 




Writing heat balance equation for tank 1 


q -q a ~9i 


A 


dh x 

dt 


9i=— q a = 

K, 


R. 


= 9 ~ 


A 


A 

A 


A 


dh x 

dt 


writing the balance equation for tank 2 


d\ di 


A 


dh 7 
dt 


A ^2 

R x R 2 


A 


dh 2 

dt 


writing steady state equations 


q s - 


R, 


h x s 

~A 


o 











h\S h 2 s _ 
i?i R 2 

writing the equation in terms of deviation variables 


Q~H l 


A 


1 1 
+ — 
R 


= A, 


i J 


dH y 

dt 


= a dH 2 _ 
R x R 2 2 dt 


taking laplace transforms 


0O)-tfiO) 


' R\+R 2 ^ 
v R iR a j 


= As H x*- 


and MAM _ MAM = Al _ S H 2 (s) 


R i R 2 


( 1 ) 

"( 2 ) 


from (1) 

MM _ 

Q(s) 


we get 


1 


A } s + 


R> + Rg 

R l R a 


H { (s) _ 

' R A " 
A +R a_ 


Q(s) ~ 

R X R aA + J 


lA+Ra 

J 






" R A ' 



if 1 ( 5 ) _ 

A +R a. 

' 7” 

_ ! AA 

Q(s) 

[v + l] 

? M 


and from (2 ) we 

get 



" R 2 R a ’ 

pO 


H x (5) _ 


l^i j 

_ T _ D J 

Q(s) 

r,5 + l](r 

2 5 + l ; 



putting the numerical values of parameters 


A (s) 

Q(s) 


_A_ 

(4 

v3 J 



8.1 A step change of magnitude 4 is introduced into a system having the 
transfer 

7 ( 5 ) _ 10 

X(5) 5 2 +1.65 + 4 

Determine (a) % overshoot 

(b) Rise time 

(c) Max value of Y(t) 

(d) Ultimate value of Y(t) 

(e) Period of Oscillation. 



I 


Given X(s)= — 7(5) 

s 


40 

s(s 2 +1.65 + 4) 


The transfer function 


is 


ns) 

X(s) 


10x0.25 

0.2(5 2 ) + ( 1 - 6 > + l) 
4 


2.5 

0.25s 2 + 0.45 + 1) 


r 2 = 0.25 ; r =0.5 

0.4 

and 2c z = 0.4 d, - ^ = 0.4 (< 1 = system is underdamped ) 

we find ultimate value of Y(t) 


405 40 

Lt 7(0 = Lt 57(5)= Lt —--= — = 10 

t ~* co s->o s->o 5(5 +1.6 + 4) 4 


thus B= 10 

now, from laplace transform tables 


7(0 - 


where 




sin(ar + </>) 



T 


,<j) =tan 



























(a) Over shoot = — = exp 
B 

thus % overshoot = 25.4 
c)thus, max value of Y(t) 

e) Period of oscillation = 


f \ 

-n% 


-;rx0.4 

U-+J 


•VO.84 


= A+B = B(0.254)+B 
= 2.54+10= 12.54 


AY 


3.427 


b) 

10 


For rise time, we need to solve 



sinCcfl 1 + (/>) 


10 fort =t r 


= e r sin(a t r + </>)= 0 

_0.4r r 

= e~ °- 5 sin(1.833t r +1.1589)= 0 
solving we get t r = 1.082 
thus 

SOLUTION: % Overshoot = 25.4 
Rise time = 1.0842 
Max Y(t) = 12.54 
U(t) Y(t) = 10 

Period of oscillation = 3.427 


= 0.254 


Comment : we see that the Oscillation period is small and the decay ratio 
also small = system is efficiently under damped. 



0, 10 ft 3 of wateris 
in level in both tanks 
at which each 


Aj = A 2 = 10 ft 3 

Ri = O.lft/cfm R 2 = 0.35ft/cfm. 



8.2 The tank system operates at steady state. At t = 
added to tank 1. Determine the maximum deviation 
from the ultimate steady state values, and the time 
maximum occurs. 


As the tanks are non interacting the transfer functions are 

//(.sp _ K, 0.1 
Q(s) r l s + 1 (s +1) 

H 2 (s) _ R 2 0.35 

Q(s ) “ (r,j +1 )(t 2 s +1) ~ (s + 1)(3.5j + 1) 

Now, an impulse of d(t) = 10 ft 3 is provided 


Q(s)=\0=H x {s) = ^—=e- t 
s +1 


and HJs) = - : -=-■=— : - 

2 (s + l)(3.5s + l) 3.5s 2 + 4.5s+ \ 


Now r 2 = 3.5= t = 1.871 


2^r = 4.5 = ^ = —=1.202 
2 r 


thus, this is an ovedamped system 





















Using fig8.5, for £ = 1.2, we see that maximum is attained at 


— = 0.95,t = 1.776min 


r 


And the maximum value is around r 2 = 0.325 Y 2 (t) = 0.174 
= H 2 (t) = 0.174x3.5 = 0.16ft 

thus max deviation is Hi will be at t = 0 = Hi = 1 ft 

max deviation is H 2 will be at t = 1.776 min = H 2 Max = 0.61 ft. 

comment : the first tank gets the impulse and hence it max deviation turns out 
to be higher than the deviations for the second tank. The second tank exhibits an 
increase response ie the deviation increases, reaches the H 2Ma x falls off to zero. 


8.3 The tank liquid level shown operates at steady state when a step 
change is made in the flow to tank l.the transaient response in critically 
damped, and it takes 1 min for level in second tank to reach 50 % of total 
change. If A 1 /A 2 = 2 ,find R 1 /R 2 . calculate x for each tank. How long does 
it take for level in first tank to reach 90% of total change? 



H x ( 5 ) _ R, 


For the first tank, transfer function , x . 

Q{S) TS + 1 


^ , , HJs) 

For the second tank —-- 


R 


Q(s) (t v s + 1)(t 2 j +1) 













H 2 {s) _ R 2 

Q(s ) t x t 2 s 2 + {t x + r 2 )5 + 1 


Q(s) = - -H 2 (s) = - -^-—— 

s s t 1 t 2 s +(r,+r 2 )i' + l 


r( parameter) = J 


= AT{U 2 


For r( parameter) = 


for ^ = l,H 2 (t) = R 2 


t 


^ + /— 

v V r i r 2 y 


given, t = 1 for r( parameter) = Jz 
H 2 (t^> qo )= i? 2 (l-(0)) =i? 2 




= ^ 


1 + 


1 




A 1 A 


V V » i 2 y ; 


2 


also 2^r = r, + r 2 


^ _ i _ r i + ^" 2 


c> /l 

-= r j? = A/? = —L = -l 

2 R 1 A, 


from I 


1-0.5 = 


ft O' 1 

1 + — e r 

V rj 


= 0.5 



8.3)tf,(s) = R ' ; Hft) = R x {\-e") 

S{T l S + 1) 


t 

0.94(7 ^oo) = Rfl-e r ‘) 


t 

0.9 R 1 = Rf l-e'^96) 


t 

e ^596 = o 1; t =1.372min 


thus 

— =0.5 

r i =r 2 = 0.596 min 
t w/i = 1.372 min 


Comment : 

Small values of t 1 ,t 2 indicate the system regains the steady state quickly. Also as R l > R 2 ,the second ta 
responds more slowly to changes than first tan k. 


8.4 Assuming the flow in the manometer to be laminar function between 
applied pressure Pi and the manometer reading h. Calculate a) steady sate 
gain ,b) r ,c) E ,. Comment on the parameters and their relation to the 
physical nature of this problem. 




Assumptions: 

Cross-sectional area =a 

Length of mercury in column = L 

Friction factor = 16/Re (laminar flow) 

Mass of mercury = m rg 

Writing a force balance on the mercury 

Mass X acceleration = pressure force - drag force - gravitational force 


( AL P) d d ^= Ap x ~4f ^— Mpgh) 


L d 2 h 8 /u dh _ p x 


— r +—T—_+ h =J-L 
g dt pgD dt pg 


At Steady state, ' pg 


Ld 2 H 8 p dH 

- z —i- 

g dt pgD dt 




\k x s 2 + k 2 s + l] H(s) =k 7t p x (s) 















#iCg) 

Pj(s) (k 1 s 2 +k 2 s + 1) 


Where 


k. k -**L. k _J_. 

’ /v 2 ~ ’ ^3 ? 

8 PS D Pg 


Thus 


Hjis) _ R 

p x (s) ( t 2 s 2 + 2<^t s + 1) 


Wliere 




2£r 


8 // 


Now 



<0# 


8 // 1 
pgD 2 r 



-1 


Steady state gain 
Lt G(s ) = i? = —; 

S^° pg 

Comment: a) T is the time period of a simple pendulum of Length L. 

b) c is inversely proportional to T , smaller the T ,the system will 
tend to move from under damped to over damped characteristics. 


8.5 Design a mercury manometer that will measure pressure of upto 2 
atm, and give responses that are slightly under damped with £ = 0.7 

Parameter to be decide upon : 


.a) Length of column of mercury 
.b) diameter of tube. 




Considering h max to be the maximum height difference to be used 


Pi = Pg h max =/ *max = 

*1^=1.5\W, 


2*1.01325*10 5 

9.81*13600 


Assuming the separation between the tubes to be 30 cm, 
We get an additional length of 0.47 m; 

Which gives us the total length L= 1.5176.47 
L = 2 M 


Now, £ =0.7 = 


4 // 

PgD 


r i —\ 




= 0.7 


4/C 


18 4*1.6*10 3 * 9 ' 81 


D = -- 

O.lpg 0.74*13600*9.81 

= 0.00015 


= 1.5*10 


As can be seen, the values yielded are not proper, with too small a diameter 
and too large a length. A smaller % value and lower measuring range of 
pressure might be better. 















8.6 verify that for a second order system subjected to a step response, 


7 ( 0 = 1 - 


Vi 7 ? 


rg ‘ sin 


Vi 7 ? 


-,i/iW 


.t 

|— + tan 

r £ 


With £ <1 


W = 


5 (r 5 + 24 r ra +1) 


r 2 ^ +2£73 , + 1 = (s-s 1 )(s-s 2 ) 


_ p )p 2 - \ 

where s l = —— + —- =-a + b 

T T 


-f . 


S 2 = ^ + 


=-a-b 


' 1 ^ 


Y(s) 

Y(s) 


\ r 2 J 


sis-s^is-s^ 


1 


A B C 

— +-+ - 


S (s-s x ) (s-s 2 ) 


^ W = -y 


A B C 

—i-+ - 


S {s-s x ) (, s-s 2 ) 


A(s 2 - (s, + S s )s + 5j5 2 ) + Bs(s -s 2 ) + Cs(s -s x ) = 1 


A+B+C =0 



- A(s j +s s )~ Bs 2 - Cs x = 1 


As y s s = 1; A = —-— = B + C = -- 


Bs 2 + Cs 2 =- 

s, 


Cs x + Cs 2 = 


5i +5 2 1 _ 1 


1 


1 1 


1 


= C= - =B = — 

^2(^2 — O) ^2(^2 — O) O'O ^1(^2 — O) 


Y(s)=- 


f , A 


1 


\ S l S 2 J 


111 11 
+ - 


8.6 7(0 = -4 


^1(^2 =o)’(s-o) ^2(^2-o)’(^-^) 

1 1 .V 1 


-e ' + 


^1^2 ^1(^2 7 ) ^2(^2 ^1) 


r ‘Si^2 


= 1 


1 


r (Jj-Jj) 


-e s '‘ - —e' sv 


J'(0 = 1 -— l — 1 

(*^2 7) 


S.t 

s 2 e ' - s,e 




F(r) = 1 4 - . r [v - V f | 



T 


7(0 = 1 - 


zje 


2 #^! 


[(-a - jb)(cosbt + jsmbt)-(-a + jb)(cosbt - j sinbt] 


Y(t) = 1- *\ e [- 2 jb(b cos bt + a sin bt )] 

2^1 


7(0 = 1 


e r 


l/T 7 ! 


^I — c ' cos(a/) + ^sin(ttf) 



(j) = tan 
verified 

8.14 From the figure in your text Y(4) for the system response is 
expressed 


viHZ 


b) verify that for B, = 1, and a step input 




Y(s) = i-5-5^ 


S T~S~ + TS + 1 


1 A Bs + C 

7 (.s')- - = — +-- 

■s(zs + l) 2 B (r + 1) 2 


A(t 2 s 1 + 2ts + 1) + Bs Cs= 1 


At 2 +5 = 0 


2 At +C = 0 


A=l; B = t 2 ;C = 2t 


Y(s) = i _ r ( n + *) + r 


s v 


(s + lf 


Y(s)A— r 


s (ts + l) (rs +1) 2 


,, & •' | 

Y(t) = 1 - e r — te 
r 


7(0 = 1 - (1 + —)fe 
r 


proved 


c) for c >1, prove that the step response 


7(0 = 1 - e T [cosh(etf) + P sinh(cKf)] 





Now Y(s) 


Hr 2 

s(s-B x )(s-s 2 ) 


Where 



T T 

from 8.6(a) 


m 


iii i i i i 

—--i- 

r 5^2 S 5j( 5 2 - 5j) (5 - Sj) s 2 (s 2 ~ Sy) (s - S 2 ) 




7(0 = 1 + e 


4 


f e*-e-*^ 




r e a ‘ +e~ a ‘^ 


Y(t) = l-e T [cosh(ctf) + j3 sinh(etf)] 

8.7 Verify that for a unit step-input 


(1) overshoot = exp 


(2) Decay ratio = exp 


vVw 1 / 

f -2 ' 


VY? 


For a unit step input the response (£<1): 


F(t)=l- 


Vi 7 ? 


- Sin 


i/Yj 


2 —+ tail 1 


£ 




J) 


(1) we have to find time t where the maxima occurs 
= dY/dt = 0 


dY) & 


dt Tyl 1 - 


- Sin 


-Cos 


Vw 


2 - + tan- 1 


£ 


= tan 


VY? 


: — + tail 1 


'Vt 


£ 


f 

i v 

Y? 

V 

£ 

: 2 ^ 

\ 

— 

= 0 

?. 

/ 

Vi" 

-# 2 


-Vt 




yy 




ve? 


= UK 


for maxima 


VY? 

£ 


= 2/; n 



t = 




8.8 Verify that for X(t) =A sin rot, for a second order system, 

Y(t ) = = ^sin (cot + </>) 


<j> = — tan 


2£,cox 
1 - (cor) 1 


Y(s) 


A _ 1 

( 5 2 + CO 2 ) ( r 2 s 2 + 2%TS + 1 ) 


r(») - # 


A 1 B' C 1 

-1-1-b 

(s-jco) s + jco (ar-Jj) 


Now as t —> co,Y(t) = A u coscot + B u sin&tf 


zV 

0-s 2 ) 


Where A n = A 1 + B l 
B n =j(A l -B 1 ) 

to detennine /l 1 , put s = jco-jco in the order 


A‘ = - J - - B' = - J - - 

2 a>(ja> -s 1 )(J(o-s 2 ) 2co(jco + s l )(jco + s 2 ) 


A 11 


j 

2 co 


_1_ 

Oh + j<o)(s 2 + j(0) 


_1_ 

(jco- Sl )(jco-s 2 ) 





AI-(cot) 2 ) 

(1 - (cot) 2 ) 2 + (I^cot) 2 



Thus tan </> 


A 11 _ - long 
B u _ 1 - (cot) 2 


And, A New = 


A 


^(W-ion) 1 ) 2 +(2con>) 2 


(using + B n = A 


New 


Thus, Y{t ) 


, Sin(cot + (/)) 

v((l - c ^^) 2 ) 2 + 


proved 


8.9) If a second- order system is over damped, it is more difficult to 

determine the parameters £, & r experimentally. One method for determining 
the parameters from a step response has been suggested by R.c Olderboung 
and H.Sartarius (The dynamics of Automatic controls,ASME,P7.8,1948),as 
described below. 

(a) Show that the unit step response for the over damped case may 
be written in the form. 


Where ri and r 2 are the roots of 
T 2 S 2 + l^TS + 1 = 0 

(b) Show that s(t) has an inflection point at 

f _ ln(r 2 /r,) 

(r 2 -r i) 


© Show that the slope of the step response at the inflection point 


d(s ) 
dt 


=s\t,) 


Wliere, .v 1 (f) =-i\e r,t = -r 2 e r -‘ 



(d) Show that the value of step response at the inflection point is 

$*(/,.) =1 + and that hence I - S ^‘^ =-—-— 

r x r 2 s (t,) r x r 2 


(e) on a typical sketch of a unit step response show distances equal to 

AO AO 

(f) Relate g & r to r x & r 2 


(a) G(s) = 


1 


t 2 s 2 + 2 %zs +1) 


f 

2 , 1 

(2£ 

1 1 1 

s + 


V 

V T J 



( s ~ r \)( s ~ r i) 


= Y{s) = 


s(s - r, )(.S' - r 2 ) 


1 ABC 

- —-1-1- 

s(s - r x )(s - r 2 ) s s-r x s-r 2 

1 = A(s - r x )(s ~r 2 ) + Bs(s - r 2 ) + cs(s - r x ) 
Put s= 0 = Arir 2 =1 ; A = t 2 

Put s = r x = Br x {r x -r 2 ) = l;B = ——^-- 

r i( r \ ~ r i) 

s = r 2 = Cr 2 (r 2 -r,) = l;C = ——--- 

r 2 {r 2 -r x ) 


f _2 


Y(s)=— 


- + - 


1 


- + - 


1 


s r x (r x - r 2 )(s - r x ) r 2 (r 2 - r x )(j - r 2 ) 


V 




(b) For inflection point, = 0 & ^ = 0 

dt dr 

ds _ r l r 1 (e r2 ‘ -e r -) 
dt r x - r 2 


d 2 s _ r l r 1 (r 2 e r2 ' -r ] e ri ‘) 
dt 2 i\ -r 2 





Also 


ds(t) 

dt 


ds{t ) 
dt 


ds(t ) _ r x r 2 (e rit - e r ' t) 

dt t=ti (r x - r 2 ) 


r,t f \ 

r i r 2 ~ g r L _ l 

d\- r 2)V2 , 


(d) 




= s(t ( ) 



Now 




proved. 

8.10 Y(0),Y(0.6),Y(oo) if 

= 1 25(, + l) 

5 (s 2 +25 + 1) 




m 


i 



Y(t) = l+5.0.3e _t sin (4.899t)-1.02e' t sin(4.899t+1.369) 


Y(0)= 1-1=0 

Y(0.6) = 1+0.561+0.515 

Y(qo) =1 

Comment : as we can see ,the system exhibits an inverse response by 
increasing from zero to more than 1 and as t tend to oo,will reach the 
steady state value of 1 . 


8.11 In the system shown the dev in flow to tank 1 is an impulse of 
magnitude 5. Ai = 1 ft 2 , A 2 = A 3 = 2 ft 2 , Ri =1 ft/cfm R 2 = 1.5 ft/cfm. 


(a) Determine Hi(s), H 2 (s), 
H 3 (s) 





















from tank 2, 


H 2 (s) 

Q(s) 


(r,s + l)(r 2 s +1) 


1.5 

(,s + l)(3s +1) 


for tank 3, 


<h -q e 


= A 3 


dh 3 

dt 


q 3 = q c {const) = Q 2 


A 


dh 3 

dt 


H, 

R 2 


a 3 


dh 3 

dt 


thus, 


A 3 SH 3 (s) = 


H 2 (s) 

r 2 


H 3 (s) _ 1 
H 2 (s ) 3s 


A 3 SH 3 (s) = 


H 2 (s) 

R 2 


HAs) 

H 2 (s) 


3s 


H 3 {s) _ 0.5 

Q{s) ~ s{s + l)(3s +1) 


8.11© HAs) =- 

1 (J + 1 ) 


//, (0 = 5e~* 


if, (3.46) = 0.155.4 


H 2 {s) _ 1.5 

Q(s) ~ 3s 2 + 4s +1 


Q{\)=5 = H 2 {s) 


1.5 

3s 2 + 4s +1 



r = V3 


2£r = 4 


£ 


_4_ 
2 T 


4 

2V3 


1.155 


from fig 8.5 
^ = \.155 and — = % = — 

T T 


3.46 


= 2 


tH 2 {t) =0.265X7.5 


H 2 (t) 


0.265X7.5 

T 


1.147 


7 / 3 (.y) _ 0.5 

0(s) ~ s(3s 2 + 4s +1) 


Q(s) = 5 = H 3 (s) 


r = V3 


2 

7J 


2.5 

s(3s 2 +4^ + 1) 


from fig 8.2 at — = 2,^ = 1.155 

T 


Y(t) =0.54 


H 3 (t) =0.54*2.5= 1.35 



8.12 sketch the response Y(t) if 7 ( 5 ) = —-- 

(5 +1.2 s +1) 


Determine Y(t) for t = 0,1,5, qo 

= e' 2 * = ^ = 1 

(s 2 +1.25 + 1) (5 + 0.6) 2 + (0.8) 2 0.8 (5 + 0.6) 2 + (0.8) 2 

7(0 = 1,25e”' 6( *“ 2) sin(0.8(t - 2)) t>2 

for t= 0 7(0) = 0 
t = 1,7(1) = 0 
t = 5,7(5) = 0.14 
t = qo , 7 (00) = 0 



Problem 8.13 The system shown is at steady state at t =0, with q = 10 
cfm 

Ai = lft 2 ,A 2 =1.25ft 2 ,R 1 = 1 ft/cfm, R 2 = 0.8 ft/cfm. 

a) If flow changes fro 10 to 11 cfm, find H 2 (s). 

b) Determine H 2 (1),H 2 (4),H 2 (qo) 

c) Determine the initial levels hi(0),h 2 (0) in the tanks. 

d) obtain an expression for H 2 (s) for unit step change. 










Writing mass balances, 

q - fe—^ =A ] d ^ (for tan k 1) 
R , dt 


At steady state q 2 -—- — = h ls ~h 2S 

R 


Also for tank 2 
(/?! - h 2 ) h 2 dh 2 
R x R 2 2 dt 

At steady state ^ ~ h ^) = !h^ = u =0.8*10=8 

1 0.8 


/?i S =18 


C) h l s = 18 ft h 2 (0) -8ft 


The equations in terms of deviation variables 


0 - 0 ,= 

01-02 


A. 


dH l 

dt 

dH , 


dt 


where 

02 = 


01 


r 2 


h x -h 2 


H 2 (s) _ R, _ 

Q(s ) t 1 t 2 s 2 +(t 1 + t 2 +A 1 R 2 )s + \ 


0.8 

s~ + 2.8^ +1 






















H 2 (s) 


0.8 

5 ( 5 2 + 2.85 + 1 ) 


(Ans 8.31(a)) 


Step response of a second order system 
r 2 = 1 = r = 1 

2£r = 2.8;£ = ^ = 1.4 


d)t = 1 = - = 1 ; H 2 (t) = 0.8(0.22) = 0.11 6 ftifrom fig) 
r 

b)t = 4 = - = 4;H 2 (t) = 0.8(0.78) = 0.624/t(/rom ,/zg) 
r 

c)t —> 00 = H 2 (t —> 00 ) = 0.8 ft 

Thus 

// 2 (1) = 0.176 ft 
H 2 (4) = 0.624 ft 
H 2 (oo) = 0.Sft 

8.13(d) we have 

=a cwco 

Gl(-S)-G 2 (- S ) = ^2( J )^2( J ) 


0 ( 5 ) - 0 2 ( 5 ) = A 1 (s)H 1 (s) + A 2 is)H 2 ( 5 ) 


Qis)-A x is)H x is) 


~~ + A 2 is)H 2 ( 5 ) 



H 2 (s) 



We have 


H 2 (s) 


t x t 2 s 2 + (r x +r 2 + A x R 2 )s +1 


*2 

Deg 


R 2 H 2 (s) _ f R 2 (Q(s) - A x sH x (s) 
Deg y (1 + t 2 s ) 

Deg { 1 Q(s) ) 


H ( .s’) 1 ^ i/r 2 (,s)^ 1 Z)eg-l-r 2 5 

20) ^4 ( Deg J A x s [ Deg 

H{ s) 1 (t x t 2 s 2 + (t x + t 2 + A x R 2 )s + 1-1- t 2 s 
Q(s) sA x [ Deg 


H{ s) _ l f t x t 2 s + t x +^ 4 , i ? 2 ) s ^ _ 
20) sA x { Deg , 

H ( s) _( R 2 +R x (1 + t 2 s) 

20) yr x T 2 s 2 +(t x +t 2 + A x R 2 )s + 1 



ews e-5. 


H,d) 



M, (O 



^ =5> f-'jtlivw u 













8.14 


_ 2 (2s + 4) 

~ 7 (4s 2 + 0.8s+ 1) 


7(s) = -- } S + 2)> - 

s (4s 2 + 0.8s +1) 


r 


F(s) = 4 




1 + - 

v s) 


1 


(4s 2 + 0.8s+ 1) 


7(s) 


8 | 4 

s(4s 2 + 0.8s + 1) (4s 2 + 0.8s+ 1) 


= (step response) + (impulse response) 


Now, t = V4 = 2; 2£t = 0.8 


E, = 0.2 


also, 


t 

T 


4 

2 


= 2 


impulse response xY(t)= 4*0.63 = 2.52 (from figure) 
step response = 8*1.15 = 9.2 (from figure) 

Y(4) = 1.26+9.2 


Y(4) =10.46 



Q 9.1. Two tank heating process shown in fig. consist of two identical, 
well stirred tank in series. A flow of heat can enter tank2. At t = 0 , the flow rate of heat 
to tanl<2 suddenly increased according to a step function to 1000 Btu/min. and the temp 
of inlet Ti drops from 60°F to 52°F according to a step function. These changes in heat 
flow and inlet temp occurs simultaneously. 

(a) Develop a block diagram that relates the outlet temp of tank2 to inlet 
temp of tankl and flow rate to tank2. 

(b) Obtain an expression for T 2 ’(s) 

(c) Determine T2(2) and T 2 (x ) 

(d) Sketch the response T 2 ’(t) Vs t. 

Initially Ti = Ti = T 2 = 60°F and q=0 

W = 250 lb/min 

Flold up volume of each tank = 5 ft 3 

Density of the fluid = 50 lb/ft 3 

Fleat Capacity = 1 Btu/lb (°F) 



Solution: 


(a) For tank 1 



Input - output = accumulation 

j'T 1 

WC(Ti - T 0 ) - WC(Ti - T 0 ) = p C V ^ .(1) 

dt 

At steady state 

WC(Tis - T 0 ) - WC(Tis - T 0 ) = 0.(2) 

(l)-(2 ) gives 


iji 

WC(Ti - Tis) - WC(Ti - T ls ) = p C V — 

dt 

rrp 1 

WTi - WTi = p V — 
dt 

Taking Laplace transform 
WTi(s) = WTi(s) + p V s Ti(s) 

—^ = —-— , where t = p V / W. 

Ti(s) 1 + ts 

From tank 2 


q + WC(T! - T 0 ) - WC(T 2 - T 0 ) = p C V —- .(3) 

dt 

At steady state 

q s + WC(Ti s - T 0 ) - WC(T 2s - T 0 ) = 0.(4) 

(3) - (4) gives 


dT i 

Q +WC(T 1 -T ls )-WC(T 2 -T 2s ) = pCV — 

dt 

dT\ 

Q + WCTj - WCT 2 = p C V-- 

dt 


Taking Laplace transform 







Q (s) + WC(T!(s) - T 2 (s)) = p C V s T 2 (s) 


T 2 (s) 


1 

I + TS 


Q(s) 

wc 


+ T\(s) 


, where x 


p V/W. 


Ti ( s ) 



Q ( s ) 


(b) x = 50*5/250 = 1 min 
WC = 250*1 =250 


Ti(s) = -8/s and Q(s) = 1000/s 

Now by using above two equations we relate T 2 and Ti as below and after taking laplace 
transform we will get T 2 (t) 


T 2 (s) = 
T 2 (s) - 


Us) = 4 


1 Q(s) , 1 

1 + ts 250 (l + zs)' 

4 8 
(1 + s ) (l + s) 

1 1 

5 (1 + J)y 


Ttis) 




2 

V f 


v- 


-8 


1 1 


v 


s + 1 


T 2 (t) = (4 + 8t)e~‘ -4 


1 

(M" 




y 


(c) T 2 ’(2) = -1.29 

T 2 (2) =T 2 ’(2) + T 2s = 60- 1.29 = 58.71 °F 
T 2 ’(oo) = -4 

T 2 (oo) = T 2 ’(qo) + T 2s = 60 - 4 = 56 °F 





























Q - 9.2. The two tank heating process shown in fig. consist of two identical, well stirred 
tanks in series. At steady state Ta = Tb = 60°F. At t = 0 , temp of each stream changes 
according to a step function 

Ta’(t) = 10 u(t) Tb’(t) = 20 u(t) 

(a) Develop a block diagram that relates T 2 ’ , the deviation in the temp of tanl<2. 
to Ta’ and Tb’. 

(b) Obtain an expression for T 2 ’(s) 

(c) Determine T2(2) 

W1 = W2 = 250 lb/min 
VI =V2 = 10 ft 3 
pf = p2 = 50 lb/ft 3 
C = 1 Btu/lb (°F) 




Solution: 

(a) For tankl 


Us) 
Ta(s ) 


- , where ii= p V / W1. 

1 + TjS 


For tanl<2 

r/7-1 

WlC(Ti - T 0 ) +W2C(Tb - T 0 ) - (W1+W2)C(T 2 - T 0 )= pCV^- .(1) 

dt 

At steady state 

WlC(Ti s - T 0 ) +W2C(Tbs - T 0 ) - (Wl+W2)C(T 2s - T 0 )= 0.(2) 

(l)-(2) 

dT o 

WITi + W2Tb - W3T 2 = p V-- 

dt 


Taking L.T 

WlTi(s) + W2Tb(s)- W3T 2 (s) = pVs T 2 (s) 





where x= p V / W3. 


T 2 (s) 


1 

1 + TS 




Tb(s) 


(b) Ti = 50*10/250 = 2 min 
x = 50*5/250 = 1 min 
W1/W3 = 1/2 = W2/W3 
Ta(s) = 10/s and Tb(s) = 0/s 

Now by using above two equations we relate Ti and Ta as below and after taking laplace 
transform we will get T 2 (t) 







T 2 (s) 

TM 

TM 

T 2 {s) 

TM 
T 2 (t ) 


Y 2 T l (s) y 2 T b (s) 


(1 + s) (l + s) 


y 2 w 


(1 + s)(l + 2s) 

(1 + ,) 

5 

10 

s( 1 + ,s')(l + 2s) 

s(l + s) 

15 + 20 5 



s( 1 + s)(l + 2s) 
(15 5 

20 ) 

1 


^ s s +1 (l + 2s)J 
15 -5e~‘ -10e"^ 


(c) T 2 ’(2) = 10.64 °F 

T 2 (2) = T 2 ’(2) + T 2s = 60 + 10.64 = 70.64 °F 


Q - 9.3. Fleat transfer equipment shown in fig. consist of tow tanks, one nested inside the 
other. Fleat is transferred by convection through the wall of inner tank. 

1. Flold up volume of each tank is 1 ft 

2. The cross sectional area for heat transfer is 1 ft 2 

3. The over all heat transfer coefficient for the flow of heat between the tanks is 10 
Btu/(hr)(ft 2 )(°F) 

4. Fleat capacity of fluid in each tank is 2 Btu/(lb)(°F) 

5. Density of each fluid is 50 lb/ft 3 

Initially the temp of feed stream to the outer tank and the contents of the outer tank are 
equal to 100 °F. Contents of inner tank are initially at 100 °F. the flow of heat to the inner 
tank (Q) changed according to a step change from 0 to 500 Btu/hr. 

(a) Obtain an expression for the laplace transform of the temperature of inner 
tank T(s). 

(b) Invert T(s) and obtain T for t= 0,5,10, oo 



10 lb/hr 


T 2 


Q 


Ti 


Solution: 

(a) For outer tank 

J r r 

WC(Ti-T 0 ) + hA (Tj -T 2 )- WC(T 2 -T 0 ) = pCV 2 ^ .(1) 

at 

At steady state 

WC(Tis - T 0 ) + hA (T ls - T 2s )- WC(T 2s - T 0 ) = 0.(2) 


(1) - (2) gives 


WCTi’ + hA (Ti’ - T 2 ’)- WCT 2 ’ 


pCV 2 


c/T 2 ' 

dt 


Substituting numerical values 

dT ' 

10 Ti’ + 10 ( Ti’ — T 2 ’) — 10 T 2 ’ = 50—— 

dt 

Taking L.T. 

Ti(s) + Ti(s) - 2T 2 (s) = 5 s T 2 (s) 

Now Ti(s) = 0, since there is no change in temp of feed stream to outer tank. Which gives 
T 2 (s)_ 1 


T x (.s’) 2 + 55 





For inner tank 


Q-hA(T 1 -T 2 ) = pCVi ^ 

dt 


Qs - hA (Tig - T 2s ) = 0 
(3) - (4) gives 


Q’-hA(TY—T 2 ’) = pCVi 


dT { ' 

dt 


Taking L.T and putting numerical values 
Q(s) - 10 Ti(s) + 10 T 2 (s) = 50 s Ti(s) 

Q(s) = 500/s and T 2 (s) = Ti(s) / (2+ 5s) 

—-10T,(.s')+ l ^ ) =50.sT,(.s') 


50 ^ 
— = T\(s) 

s 


^l(s) = 


5 s ■ 


2 +5s 
1 


2 + 55 
50(2 + 55) 


+ 1 


s(25s 2 +155 + 1) 


(3) 

(4) 


T^) = 


5 5 + 


3.82 


2(2 + 5 5 ) 

I 


50 ' 


5 + 


26.18 


50/ 


T,(s) 


100 94.71 5.29 


, _3.82 1/ _26.m/ 

Tj (t) = 100-94.71 e /so -5.29e /so 
and 


Tj (t) = 200-94.71 e 


_3.82 1/ _26.m/ 

50 - 5.29 e /so 


For t=0,5,10 and oo 
T(0)= 100 °F 
T(5)= 134.975 °F 
T(10)= 155.856 °F 
T(oo) = 200°F 





Q - 10.1. A pneumatic PI controller has an output pressure of 10 psi, 
when the set point and pen point are together. The set point and pen point are suddenly 
changed by 0.5 in (i.e. a step change in error is introduced) and the following data are 
obtained. 


Time,sec 

Psig 

O' 

10 

0 + 

8 

20 

7 

60 

5 

90 

3.5 


Determine the actual gain (psig per inch displacement) and the integral time. 

Soln: 

e(s) = -0.5/s 

for a PI controller 

Y(s)/e(s) = Kc ( 1 + xfVs) 

Y(s) = -0.5Kc ( 1/s + xfVs 2 ) 

Y(t) = -0.5Kc( 1 + xf't) 

At t = 0 + y(t) = 8 Y(t) = 8 - 10 = -2 
2=0.5Kc 
Kc = 4 psig/in 

At t=20 y(t) = 7* Y(t) = 7-10 = -3 
3 = 2(1 + tf 1 20 ) 

ti = 40 sec 

Q-10.2. a unit-step change in error is introduced into a PID controller. If Kc = 10 , n = 1 
and Id = 0.5. plot the response of the controller P(t) 

Soln: 

P(s)/e(s) = K c ( 1 + id s+ 1/ iis) 

For a step change in error 




P(s) = (10/s)(l + 0.5 s + 1/s) 
P(s) = 10/s + 5 + 10/s 2 
P(t)= 10 + 5 8(t)+ lOt 



Q - 10.3. An ideal PD controller has the transfer function 
P/e = K c (x D s + 1) 

An actual PD controller has the transfer function 
P/e = K c (x D s + 1) / (( x d /( 3) s + 1) 

Where (3 is a large constant in an industrial controller 

If a unit-step change in error is introduced into a controller having the second 
transfer function, show that 

P(t) = Kc ( 1 + A exp(-(3t/ t d )) 

Where A is a function of (3 which you are to determine. For (3 = 5 and Kc =0.5, 
plot P(t) Vs t/ Id. As show that (3 ^ oo, show that the unit step response approaches that 
for the ideal controller. 

Soln: 

P/e = Kc (x D s + 1) / (( x D /(3) s + 1) 

For a step change, e(s) = 1/s 

P(s) = Kc s( t d s + 1) / (( t d /| 3) s + 1) 






= K 


c 


1 + (J3-I)e 



So, A = (3-1 

P(t) = 0.5 ( 1 + 4 exp(-5t/ xd)) 



As (3 ^ oo then xd/(3 4 0 and 

P/e = Kc ( id s + 1) / (( xd/(3) s + 1) becomes 

P/e = Kc ( xd s + 1) that of ideal PD controller 


Q - 10.4. a PID controller is at steady state with an output pressure of a psig. The set 
point and pen point are initially together. At time t=0, the set point is moved away from 





the pen point at a rate of 0.5 in/min. the motion of the set point is in the direction of lower 
readings. If the knob settings are 
Kc= 2 psig/in of pen travel 
Ti = 1.25 min 
Id = 0.4 min 

plot output pressure Vs time 

Soln: 

Given de/dt = -0.5 in/min 
s e(s) = -0.5 

Y(s)/e(s) = K c ( 1 + id s+ 1/ xis) 

Y(s) = -( 1/s + 1/ TiS 2 + Td) 

Y(t) = -( 1 +1/1.25 + 0.4 5(t) ) 

Y(t) = y(t) - 9 = - ( 1 +1/1.25 + 0.4 8(t) ) 
y(t) = 8 - 0.8 t - 0.4 8(t) 



Q - 10.5. The input (e) to a PI controller is shown in the fig. Plot the output of the 
controller if Kc = 2 and Ti = 0.5 min 








0.5 


-0.5 


e(t) = 0.5 ( u(t) - u(t-l) - u(t-2) + u(t-3)) 


e(s) 

= (0.5/s) ( 1 - e' s 

- e' 2s + e' 3s ) 

P(s)/e(s) = K c ( 1 + (1 

/ Ti s) ) = 2 (1+ 2/s ) 

P(s) 

= ( 1 /s + 2/s 2 ) (1 

- e' s - e' 2s + e' 3s ) 

P(t) 

= 1 +2t 

0 <t< 1 


= 2 

1 <t<2 


+-» 

(N 

1 

II 

2 < t < 3 


= 0 

3 < t < oo 


Q - 12.1. Determine the transfer function Y(s)/X(s) for the block diagrams shown. 
Wxpress the results in terms of Ga, Gb and Gc 








Soln. 

(a) Balances at each node 

(1) = GaX 

(2) = (1) - Y = GaX - Y 

(3) = Gb(2) = Gb(GaX - Y) 

(4) = (3) + X = Gb(GaX - Y) + X 
Y = Gc(4) = Gc (Gb(GaX - Y) + X) 

= GaGbGcX - GbGcY + GcX 

Y _ GcjGaGb +1) 

X~ 1 + GbGc 

(b) Balances at each node 

(1) = X - (4) 

(2) = Gb(l) = Gb( X - (4)) 

(5) = GcX/Ga 

(3) = Gc(2) = GbGc( X - (4)) 


(4) = (3) + (5) .5 

= GbGc( X - (4)) + GcX/Ga 
Y = Ga(4) 

From the fifth equation 

(4) = GbGcX - GbGc(4) + GcX/Ga .6 


Q GaGbGc + Gc)X 
(1 + GbGc)Ga 

From the sixth equation 

Y _ (GaGb + \)Gc 
X~ (l + GbGc) 

Q- 12.2 

Find the transfer function y(s)/X(s) of the system shown 






Soln: 

Balance at each node 

(1) = X-Y .(a) 

(2) = (1) + (3) .(b) 

(3) = Gl(2) where G1 = l/(xis + 1) .(c) 

(4) Y = G2(3) where G2 = 0.5/(x lS /2 + 1) .(d) 

From (d) and (c) 

Y = (2)G1G2 

= G1G2 (X - Y + (3) ) .(e) 


Also from (b) and (c) 

(3) = Gl((l) + (3)) 

(3)(1 - l/(xis + 1)) = l/(ns + 1) 

(3) xis= 1 

(3) = l/(x lS ) = (X-Y) / (x!s) 

Substitute this in (e) 

Y = -- 1 + — (X - Y) 

(T lS + \)C% +1)L 

Y _ 1 

X t^s 2 +2t x s + \ 


Q - 12.3. For the control system shown determine the transfer function C(s)/R(s) 













Soln. 

Balances at each node 

(1) = R-C .(a) 

(2) =2(1) = 2(R-C) (b) 

(3) = (2) - (4) = 2(R - C) - (4) .(c) 

(4) = (3)/s = (2(R - C) - (4))/s .(d) 

(5) = (4) - C (e) 

C = 2(5) (f) 


Solving for (4) using (d) 

s (4) = 2(R - C) - (4) 

(4) = 2(R - C) / (s +1) 

Using (e) 

(6) = 2(R - C) / (s +1) - C 
C = 2 — (R-C)-C 

U+i ’ J 

4R = C((s +1) + 4 + 2(s + 1)) 
C _ 4 

R ~ 3s + 7 


Q - 12.4. Derive the transfer function Y/X for the control system shown 












Soln. 


Balance at each node 


(1) = (5) + X 

.(a) 

(2) = (1) - (4) 

.(b) 

(3) = (2)/s 

. (c) 

cn 

II 

>■ 

.(d) 

(5) = 2 (3) 

. (e) 

(4) = 25Y 

. (f) 


From (b) 

(4) = (1) - (2) 

= (1) — s (3) from (c) 

= (1)-s 2 Y from (d) 

= (5) + X - s 2 Y from (a) 
= 2 (3) + X - s 2 Y from (e) 
= 2 s Y + X - s 2 Y 

Y = (2 s Y + X - s 2 Y)/25 
X = Y( 25 - 2s + s 2 ) 

Y _ 1 

X ~ s 2 -2s + 25 


From (f) 












13.1 The set point of the control system in fig P13.1 given a step change 
of 0.1 unit. Determine 

(a) The maximum value of C and the time at which it occurs. 

(b) the offset 

(c) the period of oscillation. 

Draw a sketch of C(t) as a function of time. 



' 


; 


<5 



C 

R 


5 

(5 + l)(25 + l) 

1 + K - * * * * 5 - 

0 + l)(2s + l) 


C_ 8 
R~ 2s 2 +35 + 9 


5 

b) C( oo) = Lt . °~ 8 -= — = 0.0889 

S ^°2s 2 +35 + 9 9 


offset = 0.0111 


c) K 




1 

2V2 
























overshoot = exp 




= 0.305 


= Maximum vslue of C = 1.0305*0.0889=0.116 
Maximum value of C = 0.116 


0.116= — 
9 


1 - 


1 


V 




re sm 




—+ tan~ 

T 


* , 


< = r tan -1 — i =1.6 

Vw 7 f 


Time at which C ma x occurs = 1.6 


2m 


(c ) Period of ociullation is T = . r = 3.166 


T =3.166 


Decay ratio = (overshoot)" = 0.093 



























13.2 The control system shown in fig P 13.2 contains three-mode controller. 

(a) For the closed loop, develop formulas for the natural period of oscillation 
r and the damping factor t,\n terms of the parameters K, r^z^and r,. 

(b) Calculate E, when K is 0.5 and when K is 2. 

(c) Do £, & r approach limiting values as K increases, and if so, what are 
these values? 

(d ) Determine the offset for a unit step change in load if K is 2. 

(e ) Sktech the response curve (C vs t) for a unit-step change IN LOAD 
WHEN k is 0.5 and when K is 2. 

(f) In both cases of part (e) determine the max value of C and the time at 
which it occurs. 



1 


k 


C r,5 + l 
a) — = - , 

R k r 


1 + t d s + 


T,Sj 


1 + 


T X S + 1 


1 + T d S + 


TlSj 


C 

u 


(r,s + l) 


1 + 


+ 1 


r f 

k 

v v 


1 + t d s + 


T lSjJ 


TjS 

k 


T,T d +- 


+ 


'k + 1 

V k 


\ 

TjS + 1 

J 


C 

R 


1 + t d s + 


TiSj 


T x s +1 + k 


1 + t d s + 


TiSj 


V 







































C k(r D z jS 2 + TjS +1) 

R (kr jT D + r,r t )s 2 + (k + l)r ; s +k 







C _ TjS 

U (k + \)t,s + (r, + kr D )TjS 2 + k 

u = 1 

S 

C= -^j- 

(£ +1)^5 + (r, + k t D )TjS + k 

Lt sf(s) = — = 0 
*-*> k 



K=2 


For a unit step change in U 


C(oo)= 0 
Offset = 0 

(e) k= 0.5, =0.671 & r =2.236 


T = 




18.95 


If k = 0.5 

C _ 2s _ s 

U 5s 2 +3s + l ’ 5s 2 +3s+ 1 

If k = 2 

C _ 0.55 _ 0.5 

U ~ 2s 2 +1.55 + 1 ’ ~ 2s 2 +1.55 + 1 


In general C(t ) 


11 "7 ' /I -£ 2 " t 

—srnVW - 

ww r 


The maximum occurs at t 




2 


£ 


If k = 0.5 t max = 2.52 C max =0.42 
If k = 2 W= 1.69 C max =0.19 




13.3 The location of a load change in a control loop may affect the system 
response. In the block diagram shown in fig P 13.3, a unit step chsange 
in load enters at either location 1 and location 2. 

(a) What is the frequency of the transient response when the load enters at 
location Z? 

(b) What is the offset when the load enters at 1 & when it enters at 2? 

(c) Sketch the transient response to a step change in Ui and to a step 
change in U 2 . 


(J kr 7 ^ ^ 'Y^f c 

I 



U l = -;U 2 = 0 
s 


C S(R-C) + U ,) 


f 2 >1 

r 1 ^ 

v 2 s +1 j 

k 2s +1 j 


c 



























R = 0 


2 

C _ (2s + 1) 2 _ 2 

U ] ~ , 2 ~ As 2 + 45 + 11 

X (2s +l) 2 +1 


C _ 2 

C/j ~ % 2 +4^ + 11 





1 

VTT 


Frequency = / = — = — ^ ^ = — — = 0.2516 
T 2k r 2n 2 


C(oo) = 2/11 
Offset = 2/11 =0.182 


Ui=0;U 2 =l/s 


5x 


V 


25 + 1 


A/ 


A 


25 + 1 


C 


R=0 



1 


c _ 2s +1 
U 7 10 1 

25 + 1 25 + 1 

C _ 25 + 1 

U 2 4s 2 +45 + 11 

C(oo) = 1/11 

Offset = 1/11=0.091 


a) if U l frequency = 0.2516 

5 

if U 2 = —; frequency = 0.2516 

5 

b) if C/j =-; frequency = 0.182 

5 

if U 2 =-; frequency = 0.091 

5 












13.5A PD controller is used in a control system having first order process and a 
measurement lag as shown in Fig PI3.5. 

(a) Find the expressions for £, and x for the closed -loop response. 

(b) If Ti = 1 min, x m = 10 sec, find Kc so that '%= 0.7 for the two cases: (1) 
td =0,(2) td =3 sec, 

(c) Compare the offset and period realized for both cases, and comment on the 
advantage of adding derivative mode. 



K c (l + t D s) 

a) c_ (y + i) 

R | ^ K c (l + t d s ) 


+ l)(T m S + 1) 
























c_ = _ K c j\ + T D s){\ + T m s) _ 

R T l T m s2 + ( T l + T ,» + K C T D )•'>’ + (K c + 1 ) 





1 T l +T m +k c T D 



b) t { = lmin; r m =10 s;^ = 0.7 




Forz D - 0;C(co) = 0.76 ;offset - 0.24 
forz D = 3s; C(co) = 0.84; offset = 0.16 


for t d 


for t d 


0; period 


3 s; period 


2/rx. 

( 600 

V 4.167 

Vi 

-^ 2 

2 7Z^ 

1 600 

J 6.255 


aA"( 0-7) 2 


105.57 = period 


86.175 = period 


Comments: 


Advantage of adding derivative mode is lesser offset lesser period 

13.6The thermal system shown in figP 13.6 is controlled by PD controller. 
Data ; w = 250 lb/min; p = 62.5 lb/ft 3 ; 

Vi = 4 ft 3 ,V 2 =5 ft 3 ; V 3 =6ft 3 ; 

C = 1 Btu/(lb)(°F) 


Change of 1 psi from the controller changes the flow rate of heat of by 500 
Btu/min. the temperature of the inlet stream may vary. There is no lag in the 
measuring element. 

(a) Draw a block diagram of the control system with the appropriate transfer 
function in each block.Each transfer function should contain a numerical values of 
the parameters. 

(b) From the block diagram, determine the overall transfer function relating the 
temperature in tank 3 to a change in set point. 

(c) Find the offset for a unit steo change in inlet temperature if the controller 
gain Kc is 3psi/°F of temperature error and the derivative time is 0.5 min. 



Co^ "I ^ X 


F : 1 P' 3 *^ 

1 <3 






















WT 0 C + q = pCV l (T l -T 0 ) + WT X C 


WT ] C = pCV 2 (T 2 -T t ) + WT 2 C 


wt 2 c= pCV 3 (r 3 -t 2 )+wt 3 c 


T 0 ( WC+pCV 1 ) + q = T { (W C + pCV l ) 

T i=T 0 + --- 

1 0 WC + pCV l 

Ti=T 2 =T 3 

T = T +_-_=> T (s) = __ 

3 0 wc+pcv l 3 J (: wc+pcvjs 



13.6 (b) 


T’ ( 5 ) _ kc(l + T ° S ^ 0 +1)(1.25^ +1)(1.55 +1) 


\ + k c (\ + T D s )-—- 

C (5+ 1)(1.255 + 1X1.55 + 1) 


T,(s) __ 2k c {\ + t d s) _ 

R(s ) (5 +1)(1.8755 2 + 2.15s +1) + 2k c {\ + t d s) 

T',(s) __ 2k c (\ + T D s) _ 

R(s ) 1.8755 3 + 4.6255 2 + (3.75 + 2k c z D )s + 2k c +1 





































c)k c =3;r £) = 0.5,offset-?, t 0 ( 5 ) 


s 



r.’(s) ' 5 - >0 1.875s 3 + 4.6255 2 + (3.75 + 2k c r D )s + 2k c +1 


2k c +1 7 

Offset =0.143 

13.7 (a) For the control system shown in fig P 13.7, obtain the closed loop 
transfer function C/U. 

(b) Find the value of Kc for whgich thre closed loop response has a ^ of 


2.3. 


(c) find the offset for a unit-step change in U if K c = 4. 



U 


h'$ ¥ n 7 


v 



Js 


c 


\ 


c 


s 


U , K c s +1 

1 + ——.- 

s 0.255 +1 


























0.25s +1 


C 


U 0.25s 2 +s + i6 c (s + l) 


C 


s + 4 


U s 2 +4(K c +l)s + 4K c 


b) ^=2.3 


1 K c +1 

T= 7 - ;2^r = - 


^K~c 


K, 


c 


1 2x2.3=*‘ +1 




K r 


K r +1 
c =2.3 


yj^-c 


K c =2.952 


C) K C =4,U = 1/s 


= C =—x 


s + 4 


s s + 20s+ 16 


C(°°)= - 

4 

offset = 0.25. 


13.8 For control system shown in Fig 13.8 

(a) C(s)/R(s) 

(b) C(oo) 

(c) Offset 

(d) C(0.5) 

(e) Whether the closed loop response is oscillatory. 



U~o 



F? PI3 


4 


(a) 


C 

/? 


s(5 +1) 

1 + ^— 
.S'(.S’ + 1) 


C_ 4 
i? s 2 +s + 4 

b) C(qo) =2*1=2 

C(oo) =2 

C) offset = 0 


d) r = -;2£r = -=>£ = - 

2 4 4 

Vf5 t _i /t"t 

-h tan V15 

4 r 

14 ) 

+ tan -1 Vf5 
4 

C(0.5)=0.786 

.e) £<1, the response is oscillatory. 


= C(0.5) = 2 


1 - 


Vl5~ 


sin 


C(t) 


= 1 


sin 


A i \ 


13.9 For the control system shown in fig P13.9,determine an expression for 
C(t) 




































if a unit step change occurs in R. Sketch the response C(t) and compute 
C(2). 



C 

R 


1 + 


. + *1 

V s) 


1 + 


C_ 5 + 1 
R ~ 2s + r 



s 


c 


5 + 1 1 

-—-h 

5(25 + 1 ) 5 


-1 

25 + 1 


C(t) = l--e 
2 


t 

2 


C(2) = 0.816 


13.10 Compare the responses to a unit-step change in a set point for the 
system shown in fig P13.10 for both negative feedback and positive feedback.Do 
this for K c of 0.5 and 1.0. compare the responses by sketching C(t). 




















-ve feed back: 



+ve feed back 


C 

R 


K r x 


5 + 1 


1 ~K r 


5 + 1 


c= 




5(5 + (i-a: c )) 


For Kc = 0.5 , response of -ve feed back is 


C = 


1 


■=3+. 


2 

3 


5(25 + 3) 5 25 + 3 


1 1 1 M 


response of +v feed back is 




















1 


c= 


s(2s +1) 


1 -2 

—- 1 - 

s 2s +1 


C(t)-\-e 2 

For Kc = 1, response of -ve feed back is 
1 

s(s + 2) s s + 2 


C(t)-—-—e 
2 2 


■It 


response of +ve feed back is 



C(t) = t 


14.1 Write the characteristics equation and construct Routh array for 
the control system shown, it is stable for (1) Kc= 9.5,(ii) KC =11; (iii) Kc= 
12 


^ c 



Characteristics equation 






0 


(5 +1)0 + 2)0 + 3) 
or 0 +1)0 + 2)0 + 3) + 6Kc = 0 
( s 2 + 6s +1 Is + (6 + 6Kc) = 0 
s 3 + 6s 2 +1 Is + (6 + 6Kc) = 0 


Routh array 

s 3 1 11 

s 2 6 (6 + 6Kc ) 

s 6(1 + Kc) 


For Kc=9.5 

= 10-(Kc)= 10-9.5=0.5>0 therefore stable. 
For Kc=ll 

= 10-(Kc)= 10-11=-1<0 therefore unstable 
For Kc=121 

= 10-(Kc)= 10-12=—2<0 therefore unstable 


14.2 By means of the routh test, determine the stability of the system shown 
when Kc = 2. 


Characteristic equation 


r 


1 + 2 


A 


1 + - 

V s) 


10 


2s 2 + 4s+ 10 


0 



(2s 2 + 45 + 10)5 + 2(s + 3).2.10 = 0 
2s 3 +4r + 105 + 405 + 120 = 0 
25 3 +4s 2 + 50s+ 120 = 0 
s + 2 s 2 + 25^ + 60 = 0 


Routh Array 

1 25 

2 60 
- 10/2 

The system is unstable at Kc = 2. 

14.4 Prove that if one or more of the co-efficient (a 0 ,ai,....a n ) of the 
characteristic equation are negative or zero, then there is necessarily an 
unstable root 

Characteristic equation : 

UqX + u^x +.+ a n — 0 

a 0 (x n + aja 0 x n ~ l +. a n /a 0 ) = 0 

a 0 (x - cq )(x - a 2 ). (x - a n ) = 0 

We have a x ,a 2 . a n < 0 


As we know the second co-efficient aj/ao is sum of all the roots 


- L = (-l ) 5 

a n 


/=1 7=1 


a , 


12 


Therefore sum of all possible products of two roots will happen twice as 
a { a 2 dividing the total by 2. 

a { aj > 0 (a i <0 a j < 0) 

And 

:.a 2 / a 0 > 0 => a 2 >0 


Similarly 







/ 

— = (-1 y ( sum of aoll possible products of j roots ) 


if j = even (-1) ' is 1 and the sum is >0 so a - / a 0 > 0 


if j = odd (-1)' is (-1) and the sum is <0 so — is again > 0 


in both case a . / a {] >0 
so a j > 0(forj-\, . n ) 


14.5 Prove that the converse statement of the problem 14.4 that an 
unstable root implies that one or more co-efficient will be negative or 
zero is untrue for all co-efficient ,n>2. 

Let the converse be true, always .Never if we give a counter example we 
can contradict. 

Routh array 

s + .s 1 + 2s + 3 
s 3 1 2 
s 2 1 3 
s -1 0 

s° 

System is unstable even when all the coefficient are greater than 0; hence a 
contradiction, 

14.6 Deduce an expression for Routh criterion that will detect the 
Presence of roots with real parts greater than a for any rectified <r>0 

Characteristic equation 

a 0 x" + a x x n ~ l + .+ a n - 0 


Routh criteria detennines if for any root, real part > 0. 
Now if we replace x by X such that 





.x + a =X 


Characteristic equation becomes 

a 0 (X-a)" +a 1 (X-cr)”' 1 +. + a n = 0 

Hence if we apply Routh criteria, 

We will actually be looking for roots with real part >cr rather than >0 

ci 0 x" + a x x"~ x +a 2 x"~ 2 . + a n =0 

Routh criterion detects if any root ads greater than zero. 

Is there any x = a x ,a 2 , . ,a J} . a n > 0-(1) 

Now we want to detect any root 

a j > -a 
a j> 0 
From( 1) 








x = a x ,a 2 , . aj, .a„,> 0 

implies is there any 
x = a x > 0 
x = a 2 > 0 


x - a j> 0 


x = a n > 0 

add a on both sides 
is there any 
x + <j = a l +<7>0 
x + <j = a 2 +(r> 0 


.x + cr = a j +cT> 0 


x + a = a n +(T> 0 
so. Let X = x +<7 

and apply Routh criteria to det ect any a , + cr > 0 
or a ; > -a 

14.7 Show that any complex no Si satisfying 15| < 1, yields a value of 

1 -I- e 

Z = —^that satisfies Re(Z)>0, 
l-s 

Let S=x+iy, ^ x 2 + y 2 < 1 






z = 


1 + 5 
1-5 



(l + x) + zy (1 -x) + iy 
(1 - x) - iy (1 - x) + iy 
(1-x 2 + (1 = x +1 -x)iy - v 2 ) 

1 + x 2 — 2x + y 2 

= _ l-(x 2 +y 2 ) + 2/y 

1 — 2x + (x 2 + y 2 ) 


Re(Z) 


l-2x + (x 2 +y 2 ) 


l-(x 2 +y 2 ) 


if Re(z) > 0 then 1 - (x 2 + y ) > 0 and 


1 -2x + (x 2 + y 2 ) > 0 

we have^jx 2 + y 2 < 1 
x 2 +y 2 < 1 


Ranges are -1 < x < 1 

— 1 < y < 1 Points in the unit circle 

1 - (x 2 + y 2 ) > 0 is true therefore x 1 + y 2 <1 


Now 

l + (x 2 + y 2 )-2x 

if x = -1& y = 0 then it is 4 
if x = 1 & y = 0 then it is 0 

0 < (1 + (x 2 + y 2 ) - 2x) < 4 

Re(z) > 0 


example. 

if s = (0.5 + iO.5) the system is unstable due to the real part 


s -(0.5 + z‘0.5) 



1 


= e (Cos 0.5 1 + Sin0.5t ) 

s-(0.5 + /'0.5) 

14.8 For the output C to be stable, we analyze the characteristic 
equation of the system 


1 + 


1 


1 


TjS (t y s + l)(r 2 s +1) 


x (r 3 s +1) = 0 


t i s(t 1 s 2 3 t 2 + T x s + t 2 s +1) + t 3 s + 1 = 0 


+ T I (r l + T 2 )s 1 + (Tj + T 3 )s + 1 = 0 


Routh Array 

3 

S T jT2 T j + T 3 

S 2 Ti(t 1 +T 2 ) 1 

s a 0 

5 ° 1 

= r /( r i + r 2 )( r / + r 3 ) ~ r / r i r 2 

Ti(T\+T 2 ) 


Now 

(1) t i t 2 t 1 > 0 

Since r, & r 2 are process time constant they are definitely +ve 

Tj > 0; t 2 > 0 

(2) r / (r 1 +r 2 )>0 

(3) a > 0 t i (t 1 + t 2 )(Tj + t 3 ) > TjT x t 2 



T\Ti + T l T 3 + T 2 T i + T 2 T 3 T x T 2 > 0 
r /( r i + Z" 2 ) > T x T 2 ~ T i {z l + T 2 ) 


tI >-— - r 3 

r i + r 2 

also Tj > 0 

14.9 In the control system shown in fig find the value of Kc for which 
the system is on the verge of the instability. The controller is replaced by a 
PD controller, for which the transfer function is Kc(l+s). if Kc =10, 
determine the range for which the system is stable. 



<0 A « . . /% V— 1 \0 1 A 

Characteristics equation 


1 +-= 0 

(s + 1)0 + 2)0 + 3) 

or (5 +1)0 + 2)0 + 3) + 6Kc - 0 
O 2 + 3s + 2)0 + 3) + 6 Kc = 0 
s 3 + 6 s 2 +1 Is + (6 + 6 Kc) = 0 


Routh array 


s 3 1 3 
s 2 3 1 + Kc 


s 


3- 


r 


1 + Kc 


) 




\ 


3 







For verge of instability 3 = 

Kc = 8 



Characteristics equation 


j + 10(1 +fay) _ Q 

(,y + 1)3 

s 3 +3s 2 +5(3 + 1 OKcs) + 11 = 0 

Routh Array 

5 3 1 3 + 10r D 

5 2 3 11 


s 


3(3 +10r s ) > 11 for vege 
30r s > 2 

t d > 2/30 

14.10 (a) Write the characteristics equation for the central system shown 
(b) Use the routh criteria to determine if the system is stable for Kc=4 
© Determine the ultimate value of Kc for which the system is unstable 






(a) characteristics equation 






f 5+2 l 

f 1 ^ 

| 1 

l 3 J 

125 = lj 

1(5 + 1) 


(5 2 + 5)(25 +1) + kc(s + 2) = 0 

2 s' + 3 5 2 + (1 + kc)s + 2 kc = 0 
s 3 + 3s 2 + 35 + (1 + kc) = 0 

Kc=4Routh array 

5 3 2 5 

5 2 3 8 

5 -1/3 

not stable 

3(1 + kc) - 4 kc _ 

3 

3-Kc = 0;Kc = 3 


For verge of instability 

14.11 for the control shown, the characteristics equation is 

5 4 +45 3 +65 2 +4j + (1 + A-) = 0 

(a) determine value of k above which the system is unstable. 

(b) Determine the value of k for which the two of the roots are on 
the imaginary axis, and determine the values of these imaginary roots and 
remaining roots are real. 


5 4 + 45 3 + 65 2 + 45 + (1 + k) = 0 



s 4 1 6(1 + *) 

s 3 4 4 
s 2 5 l + k 

s 4- —(1 + &) 
5 

1 l + k 


For the system to be unstable 


'l + k 

v 5 


<0 


1 < 


1 + k 
5 


k> 4 
l + k<0 
k < -1 
k> — 1 


The system is stable at -l<k<4 
(b) For two imaginary roots 

4 = ^(l + k);k = 4 


Value of complex roots 
5s 2 +5 = 0 
s = ±i 


s 2 +l s 4 + 4s 3 + 6s 2 + 4s + 5 s 2 +4^ + 5 



s 4 +0 + s 2 


4 s 3 +5s 2 
4s 3 + 0 + 45 
55 2 +5 
5s 2 +5 
0 


SOLUTION: 


4± V16-20 
2 



2 


PART 2 


LIST OF USEFUL BOOKS FOR PROCESS 
CONTROL 


1. PROCESS CONTROL BY R.P VYAS, CENTRAL TECHNO 
PUBLICATIONS, INDIA ( WIDE VARIETY OF SOLVED PROBLEMS ARE 
A VAILABLE IN THIS BOOK) 

2. ADVANCED CONTROL ENGINEERING BY RONALD.S .BURNS , 
BUTTERWORTH AND HIENEMANN. 

3. PROCESS MODELLING SIMULATION AND CONTROL FOR 
CHEMICAL ENGINEERS, WILLIAM.L.LUYBEN, MCGRAW Hil l . 

4. A MATHEMATICAL INTRODUCTION TO CONTROL THEORY BY 
SCHOLOMO ENGELBERG, IMPERIAL COLLEGE PRESS 






LIST OF USEFUL WEBSITES 


www.msubbu.com FOR BLOCK DIAGRAM REDUCTION AND 
OTHER CHEMICAL ENGG. LEARNING RESOURCES 


Readings,Recitations,Assignments,Exams,StudyMateri 
als,Discussion Group,Video Lectures 
now study whatever u want with respect to chemical 

engg. 

http://ocw.mit.edu/OcwWeb/index.htm 
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